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ABSTRACT 

The Milky Way's dwarf spheroidal satellites include the nearest, smallest and least 
luminous galaxies known. They also exhibit the largest discrepancies between dy- 
namical and luminous masses. This article reviews the development of empirical con- 
straints on the structure and kinematics of dSph stellar populations and discusses how 
this phenomenology translates into constraints on the amount and distribution of dark 
matter within dSphs. Some implications for cosmology and the particle nature of dark 
matter are discussed, and some topics/questions for future study are identified. 



1. Introduction 

Physics assigns to gravity the responsibility of forming structure on scales ranging from ter- 
restrial to cosmological. An apparent threshold arises somewhere between scales characteristic of 
star clusters and those characteristic of galaxies. The internal dynamics of gravitationally bound 
structures smaller than a few parsecs, of which globular clusters are the largest examples, are 
reasonably well described in terms of standard gravity (i.e., Einstein's general theory and/or its 
Newtonian approximation) sourced by known substances. The internal dynamics of gravitation- 
ally bound structures larger than a few tens of parsecs, of which dwarf spheroidal (dSph) galaxies 
are the smallest examples, are not. 

The ubiquity of dark matter on galactic and larger scales signifies new physics. Either there ex- 
ists an otherwise unknown substance that contributes to dynamical mass but not to baryonic mass, 
or the standard dynamical framework requires modification (or both). The 'substance' hypothesis 
is not falsifiable, but in principle it can be confirmed with the detection of non-gravitational in- 
teractions involving dark matter particles. In any case, dSphs provide the most extreme examples 
of dark matter phenomenology, with dynamical mass-to-light ratios M/Ly > 10 [M/Lv]q even at 
their centers. This fact has made dSphs the focus of intense scrutiny in the effort to understand the 
nature of dark matter. 

This article reviews the development of empirical constraints on the amount and distribution 
of dark matter within the Milky Way's dSph satellites. These results follow from the application of 
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a rich variety of analyses applied to observations conducted by many individuals and groups. All 
analyses described here are formulated within the Newtonian dynamical framework. The reader 
is welcome to interpret 'dark matter' in terms of the substance hypothesis or more generally as 
a quantification of the discrepancy between dynamical and baryonic mass. The focus here is on 
the relationship between data and constraint, and one hopes that this information translates mean- 
ingfully into alternative dynamical frameworks (as formulated, e.g. by |Bergmann|| 1968 Milgrom 
19831 |Bekenstein|2004[ pbffat|2006D . 



2. Observations 

Bright red giant stars are detectable as point sources with my < 21 mag out to distances of 
~ 0.5 Mpc. Within this range, the Milky Way's dSph galactic satellites appear as localised over- 
densities of individually resolved stars. Empirical information about the number, stellar structure 
and internal kinematics of these systems has steadily accumulated for the past eight decades. 



2.1. Census 



Shapleyl (1938a) discovered the Sculptor dSph upon visual examination of a photographic 
plate exposed for three hours with the 24-inch Bruce telescope at (what was then) Harvard's Boy- 
den observatory in South Africa. In this original dSph discovery paper, Shapley notes that Sculptor 
is visible — in hindsight — as a faint patch of light on a plate taken for the purpose of site testing, in 
1908, during a series of five exposures totalling nearly 24 hours with a 1-inch refracting telescope. 



Reporting shortly thereafter the similar discovery of Fornax, | Shapley| ( |1938b[ ) notes that the 
new type of 'Sculptor-Fornax' cluster shares properties with both globular clusters and elliptical 
galaxies, then speculates that 'At the distance of the Andromeda system these objects would, in 
fact, have long escaped discovery. There may be several others in the Local Group of galaxies; 
such objects may be of frequent occurence in intergalactic space and of much significance both in 
the census and the genealogy of sidereal systems.' 

Figure [T] demonstrates Shapley 's prescience, plotting the cumulative number and luminosities 



of the Milky Way's known dSph satellites against date of discovery publication. Harrington & 
Wilson] ( [1950| ) and [Wilson] ( | 195 5 1 ) found the next four 'Sculptor-type' (as they were then called) 
systems — Leo I, Leo II, Draco and Ursa Minor — on photographic plates taken for the Palomar 



Observatory Sky Survey with the 48-inch Schmidt telescope. Cannon et al. ( 1977 1 spotted Carina 



by eye on a plate taken with the 1 .2-meter UK-Schmidt Telescope. Irwin et al. ( |1990| ) used the Au- 
tomated Photographic Measuring (APM) facility at the University of Cambridge to detect Sextans, 
again on a UK-Schmidt plate, thereby completing the census of the Milky Way's eight so-called 
'classical' dSphs. 
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Fig. 1. — Cumulative number (top) and luminosity (bottom) of known dwarf spheroidal satellites 
of the Milky Way, versus year of publication of discovery paper ( |Shapley||1938a|b[ [Harrington & 
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The next discovery was unique in that it came via spectroscopy rather than imaging. During 
a stellar kinematic survey of the Galactic bulge, Ibata et al. (1994]) noticed Sagittarius (Sgr) in 
color-velocity space as an excess of stars following a narrow velocity distribution offset from that 
of bulge members. Subsequent observations showed that Sgr, located just ~ 15 kpc behind the 
Galactic center, spews stars along tidal tails that wrap spectacularly around the entire sky (e.g., 
Mateo et al.|[T996t |Ibata et al.|[T997t [Majewski et"aL|[2003t [Belokurov eraL][2006t [Koposov et aL 
2011a[). The tidal disruption of Sgr provides a useful tracer of the Galactic potential (e.g., [Johnston 



et al. 2005 , [Law & Majewski 2010, Penarrubia et al. 2010), but disqualifies Sgr from a simple 
equilibrium-dynamical analysis. For this reason Sgr shall not be considered further here. 

In a span of seven years, deep multi-color photometry from the Sloan Digital Sky Survey 
(SDSS) has now tripled the number of known Milky Way satellites ( [Willman et al.|2005b|a^|Zucker 



et al.||2006b|a, IBelokurov et al.||2006 


Belokurov et al.||2007 


Irwin etal.pOOV; 


Walsh etal.||2007; 


Belokurov et al. 


2008 


Belokurov et al. 


2009 


Grillmair 


2009 


Watkins et al. 


2009 


Belokurov 



et al.||2010[ ). The 17 satellites discovered with SDSS data have extended the floor of the observed 



galaxy luminosity function from My > -8 to My > -2 (bottom panel of Figure [T]), such that some 
galaxies are less luminous than some stars! Unlike their brighter 'classical' siblings, the 'ultrafaint' 
satellite^discovered with SDSS data are not apparent to the eye, even in deep images. Rather, they 
are detected only by correlating spatial overdensities with overdensities in color-magnitude space 
(e.g., [Belokurov et al.|2006| [Walsh et al.|2009| ). In order to confirm the faintest satellites, the SDSS 
catalog must be supplemented by deeper, follow-up photometry as well as spectroscopy (e.g., 
Belokurov et al.[[2009,|2010j ). Given the rate of false positives expected for candidates remaining 
in SDSS data (> 4/5 in the author's experience), the expense of follow-up observations can be 
prohibitive. However, the next generation of sky surveys (Pan-STARRS, SkyMapper, DES, Gaia, 
LSST, etc.) will almost certainly bring a new flurry of discoveries, particularly in the relatively 
unexplored southern sky. 

Although they are not considered further below, it is worth noting that Shapley was correct 
regarding the dSph satellites of M31. van den Bergh ( |1972 ) discovered the first example with 
the Palomar Schmidt telescope, using plates more sensitive than those used in the original Palo- 
mar Survey. The current census includes 27 known dSph satellites of M31. Two-thirds of this 



number were discovered in the past seven years (Zucker et al. 


2004; 


Martin et al. 2006 I 


bata et 


al. 2007 ; Zucker et al. 2007[ Majewski et al. 2007 twin et al. 


2008 


McConnachie et al. 


2008; 


Martin et al.[[2009 Richardson et al.|201 1 ; Slater et al.[[201 1 Bell et al.|201 ID, with SDSS data as 



well as photometry from the PAndAS survey conducted with the Canada-France-Hawaii Telescope 
dMcConnachie et al.|2009l). 



'As Figure [ij and the terms themselves suggest, the distinction between 'classical' and 'ultrafaint' dSphs involves 
a mixture of intrinsic luminosity with sequence of discovery. Here this distinction (which is meaningless in the 
sense that members of both classes trace smooth scaling relationships involving luminosity, size, metallicity and 
stellar kinematics) is preserved only because observational studies of these objects — for both practical and accidental 
reasons — tend to be separable along the same lines. Following common practice, dSphs known before SDSS (Carina, 
Draco, Fornax, Leo I, Leo II, Sculptor, Sextans, Ursa Minor) are referred to as 'classical', and the rest as 'ultrafaint'. 
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2.2. Stellar Structure 

Galactic dynamics is concerned with the relationship between gravitational potential and the 
distribution of stars in phase space. Observers who study dark matter in dSph galaxies must gather 
information about the positions and velocities of dSph stars. The largest dSphs subtend solid 
angles of several square degrees, making it difficult to study their stellar structure at large radius. 
Complete homogenous surveys are rare and valuable. 



2.2.1. 'Classical ' dSphs 



Hodge| ( [T961b|aH1962H1963[|1964a|b| ) used photographic plates obtained at the Palomar (48- 
inch Schmidt, 100-inch and 200-inch), Lick (120-inch) and Boyden Station (24-inch Schmidt) 
in order to study luminous structure of the six Milky Way dSphs known at the time (Sculptor, 
Fornax, Leo I, Leo H, Draco and Ursa Minor). Hodge counted stars within squares of regular 
grids overlaid on each plate, and provides this description of analog data reduction: 'Each plate 
was counted at one sitting so that uniformity would be maintained. The plates were counted to 
the limiting magnitude and each was counted once. From experience ... it was decided that the 
reproducibility of counts on a particular plate is greater than from plate to plate, so that it is better 
to count many plates each once than one plate many times' ( |Hodge|1961b| ). 

Figure [2] (top two rows of panels) displays isopleth maps that Hodge drew by connecting 
squares containing equal numbers of stars. These maps reveal that the internal structure of dSphs is 
smooth. Hodge reasoned that a well-mixed dSph requires dynamical relaxation by a process other 
than stellar encounters, as the low surface densities of dSphs imply internal relaxation timescales 
of > 10^ Hubble times. Hodge (1966) and Hodge & Michie (1969) would later suggest that en- 



counters between stellar groups might enable exchanges of orbital energy over shorter timescales 



during dSph formation. Eventually Lynden-Bell (1967) would show that the time- varying gravi- 
tational potential of a young galaxy effectively shuffles the orbital energies of its stars, generating 
'violent' relaxation without stellar encounters. More recently, [Mayer et aL (2001b) have proposed 
a mechanism specific to dSphs, demonstrating with N-body/hydrodynamical simulations that re- 
peated tidal encounters with the Milky Way can effectively transform a rotating dSph progenitor 
into a pressure-supported spheroid in less than a Hubble time. 

Hodge discovered two other structural features common to the classical dSphs. First, most 
exhibit flattened morphology, with typical ellipticities of e = 1 -b/a ~ 0.3, where a and b are semi- 
major and semi-minor axes, respectively. Second, dSph stellar density profiles decline more steeply 
at large radius than do the profiles of giant elliptical galaxies. Whereas the latter are commonly fit 
by formulae with relative ly shallow outer pro files, e.g., S(7?) = S(0)/(1 +R/df- ( Hubble] 1930 ) or 
S(i?) = S(0)exp[-fci?^/'^] ( de Vaucouleurs 1 1 948 j ) , Hodge found that classical dSphs all have steeper 



Fig. 2. — Stellar isodensity maps for the Milky Way's eight 'classical' dSphs. Top two rows: from 
Hodge's star count studies ( Hodge|1961b|a[ 1962[ 1963 , 1964a|b[ reproduced by permission of the 
American Astronomical Society). Bottom two rows: reproduced from Structural Parameters for 
the Galactic Dwarf Spheroidals, by M. Irwin & D. Hatzidimitriou, MNRAS, 277, 1354, 1995 (by 
permission of John Wiley & Sons Ltd.). 
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outer profiles that are better fit with the formula of|King|([1962[): 



E(i?) = k 



(1) 



where Rc is a 'core' radius and R^ is a maximum, or 'limiting' radius that one might expect to 



result from tidal truncation (Section 2.2.3 



Three decades later, Irwin & Hatzidimitriou (1995] 'IH95' hereafter) used the APM facility at 



Cambridge to count stars automatically on photographic plates from the Palomar and UK Schmidt 
telescopes. While confirming Hodge's findings, IH95 produced significantly deeper maps and were 
able to include Carina and Sextans, the two Milky Way dSphs discovered in the interim (Figure [2]). 
IH95 used these maps to measure the centroid, ellipticity and orientation of each dSph, and then to 
tabulate stellar density as a function of distance along the semi-major axis. While the homogeneous 
analysis of IH95 continues to provide a valuable resource particularly for comparing dSphs in the 
context of scaling relations (Section [5]), deeper photometric data sets now exist for most of the 



classical dSphs (e.g.. 


Stetson et al.||1998 , |Majewski et al.||2000 


Saviane et al.||2000 


Odenkirchen 


et al.pOOlt Palma et al.||2003 , 


Walcher et al.||2003t Lee et al.||2003 ,|Tolstoy et al.||2004; 


Coleman 


et al.|2005a|bHBattagfia et al.|2006 Westfall et al.|2006). 



In principle the structure of dSph stellar components carries information about the mecha- 
nisms that drive dSph formation and evolution. While incompatible with shallow Hubble and de 
Vaucouleurs profiles, the available data often do not distinguish the King profile (Equation[T]) from 
other commonly adopted fitting formulae — e.g., exponential and jPlummer ( |1911| ) profiles. 



S(i?) = S(0)exp[-7?/i?,] 



and 



S(i?): 



S(0) 



[\ + {R/Rn?] 



2' 



(2) 



(3) 



respectively, where i?/, ~ 1.68i?e is the projected halflight radius (i.e., the radius of the circle en- 
closing half the stars as viewed in projection). Figure [3] displays IH95's fits of King (1966) and 
exponential surface brightness profiles to the Milky Way's classical dSph satellites. 



2.2.2. Ultrafaint dSphs 



The SDSS catalog enables homogeneous studies of the structural properties of ultrafaint 
dSphs. For example. Figure |4] displays isopleth maps constructed by Belokurov et al. (2007) for 



^In fact /?K is usually referred to as a 'tidal' radius and denoted r,. The adopted nomenclature and notation avoid 
confusion with the tidal radius defined in Equation |4| 
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Fig. 3. — Stellar density profiles for the Milky Way's 'classical' dSphs. Overlaid are best-fitting 



King ( 1966, solid) and exponential (dashed) surface brightness profiles. Reproduced from Struc- 



tural Parameters for the Galactic Dwarf Spheroidals, by M. Irwin & D. Hatzidimitriou, MNRAS, 
277, 1354, 1995 (by permission of John Wiley & Sons Ltd.). 
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Coma Berenices, Canes Venatici II, Segue 1 and Hercules using both SDSS and deeper follow- 
up data. Whereas Hodge and IH95 estimated structural parameters for the classical dSphs after 
binning their star-count data and subtracting estimates of foreground densities, the low surface 
brightnesses of 'ultrafaint' dSphs are conducive neither to binning nor to foreground subtraction. 
Fortunately, neither procedure is necessary; indeed Kleyna et al.| ( |199 8) estimate structural param- 
eters for the classical dSph Ursa Minor using a likelihood function that specifies the probability 
associated with each individual stellar position in terms of a parametric surface brightness profile 
plus constant foreground. 



Martin et al. (2008) design a similar maximum-likelihood analysis that operates directly on 



unbinned SDSS data in order to estimate the centroid, halflight radius, luminosity, ellipticity and 
orientation of each ultrafaint satellite. Figure [5] demonstrates the efficacy of their method. The 
left panel maps individual stars from the SDSS catalog that are near the line of sight to Bootes I 
and satisfy color/magnitude criteria designed to select red giants at the distance of Bootes I. Panels 
on the right-hand side show maximum-likelihood estimates of each free parameter. Marginalised 
error distributions include the effects of sampling errors and parameter covariances, and can be 
used directly in subsequent kinematic/dynamical analyses (Section]?]). 

Martin et al.| ( |2008| ) show that their method recovers robust estimates even when the SDSS 



sample includes as few as tens of satellite members. Subsequent tests with synthetic data by Munoz 
et al.| ( 20TT] ) provide reason for caution, suggesting that for objects with low surface brightness, 
insufficient contrast between members and foreground can generate biased estimates of structural 
parameters. It is reassuring that deeper observations with CFHT ( [Munoz et al.||2010[ ) and Subaru 
( Okamoto et al.|20T? ) yield structural parameters for subsets of ultrafaint satellites that agree well 
with the SDSS-derived estimates of |Martin et al.H2008D . 



2.2.3. Extended Structure 



The outer stellar structure of a given dSph is determined by some combination of formation 
processes and subsequent evolution within the external potential of the Milky Way. In his structural 
analyses of the six dSphs known in the 1960s, Hodge compared the observed limiting radii, 
(Equation [TJ, to simple estimates of 'tidal' radii, r^, beyond which stars escape into the external 
potential of the Milky Way ( |von Hoerner|1957^[King|1966[ ): 



Rt 



Mr 



1 1/3 



(3 + e)MMw 



(4) 



Here, Mmw and Md are the Milky Way and dSph (point) masses, respectively, R^ is the pericentric 
distance of the dSph's orbit and e is orbital eccentricity. Assuming circular orbits and considering 



only luminous masses, Hodge (1966) noticed that while the two radii are similar for the nearest 
dSphs (Draco, Sculptor, Ursa Minor), for the three most distant dSphs (Leo I, Leo II, Fornax), 
7?K ^ 0.5r,. Hodge ( 1966) concluded that tidal forces play a significant role in shaping the outer 
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Fig. 4. — Stellar idodensity maps for four of the Milky Way's 'ultrafaint' dSphs, from the discovery 
paper of Belokurov et al.| ( 2007[ reproduced by permission of the American Astronomical Society). 
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Fig. 5. — Measurement of structural parameters for the Bootes I dSph from SDSS data (Martin 



et al. 2008 1. Left: sky positions of red giant candidates selected from the SDSS catalog. Right: 
constraints on structural parameters, from the maximum- likelihood analysis of Martin et al. ( 2008[ 
reproduced by permission of the American Astronomical Society). 
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structures of the nearest dSphs. 



Low surface brightness in the outer regions of dSphs makes extended structural studies diffi- 



cult to conduct and interpret. For example, [Coleman et al.| ( |2005a| ) use deep wide-field photometry 
to conclude that Sculptor's surface brightness profile is well described as the superposition of two 
equilibrium stellar components with different scale radii. Using an independent data set, jWestfall 
et al. ( |2006| ) fit a single profile and identify a 'break' in Sculptor's surface brightness profile, which 
makes an apparent transition from a King profile in the inner parts to a shallower power law in the 



outer parts. Westfall et al. (2006) cite this transition not as evidence for a second component, but 
rather as a signal that Sculptor is losing stars to tidal disruption. 

In many cases, the use of narrow-band filters that are sensitive to stellar surface gravity (e.g., 
DD051) can help to distinguish dSph red giants from foreground dwarf stars, providing a valuable 
boost in contrast (e.g., |Majewski et al.|2000[[2005||Parma et al.|2003||Westfall et al.'2006). Having 
used narrow-band photometry to select spectroscopic targets at large dSph radi, ,Munoz et al.| ( ,2005| ) 
and Munoz et al.| ( |2006l ) present velocities that confirm the membership of stars in Ursa Minor and 
Carina (Figure [6]) out to radii ~ 3 and ~ 5 times larger, respectively, than estimates of R^- These 
detections imply that tides can play a significant role in shaping the outer parts of dSphs (Section 



2.2.4. Structural Peculiarities of Individual dSphs 



Most kinematic analyses of dSphs proceed from the assumption that dSphs host a single, 
spherically symmetric stellar component in dynamic equilibrium (Sectio n |4]). However , real dSphs 
are all flattened <e=l-b/a< 0.7, [Irwin & Hatzidimitriou|[T995| [Martin et al.[[2008[ [Sand 
etal.|20ir[ ), and it is not clear how severely this violation of spherical symmetry affects conclusions 
regarding dSph dynamics. In fact most dSphs exhibit individual peculiarities that further violate 
the simplistic assumptions (Section[4]) employed in kinematic analyses. 

For example. Sculptor, Fornax and Sextans all display evidence for chemo-dynamically in- 
dependent stellar ^wZ^-populations ( [Tolstoy et al.[[2004^[Battagfia et al!][2006[ [Battaglia et al.|[2011 



respectively). In all three cases, a relatively metal-rich, kinematically cold population has smaller 
scale radius than does a metal-poor, kinematically hot population (Section 4.2. 2[ Figure 10). For- 



nax also has irregular stellar structure in the form of a ere scent- shaped feature near its center ( [Stet- 
son et al.[[1998[ and Figure [7} left panel), two shell-like features ( [Coleman et al.[[2004[ ) and lobes 
along its morphological minor axis ( [Coleman et al.[[2005bl and Figure [7} right panel). These fea- 
tures suggest that Fornax may have undergone a recent merger, an event that might be related to the 
presence of a young (age ~ 100 Myr), centrally concentrated main sequence in Fornax (Battaglia 



et al. 2006). Ursa Minor exhibits clumpy stellar substructure, most dramatically in the form of a 
secondary peak in its luminosity distribution, offset by ~ 20' from the central peak (Olszewski & 
Aaronson|1985 ). The region near the secondary peak is kinematically colder than the rest of Ursa 
Minor ( [Kleyna et al.|2003 ), and so may represent a bound star cluster (Section 4.2.1 ). 
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Fig. 6. — Extended stellar structure of the Carina dSph, from narrow -band photometry and spec- 
troscopy of Munoz et al. (2006} reproduced by permission of the American Astronomical So- 
ciety). The top/bottom panels show sky positions of red giant candidates confirmed as Carina 
member s/nonmembers. Ellipses mark Carina's limiting radius, (Equation [T]), determined from 
smooth fits to star count data. The extension of faint stellar structure at 7? > 37?k suggests tidal 
interaction (Section 13. 21). Munoz et al. I (12005 ) report similar results for Ursa Minor. 
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Fig. 7. — Structural irregularities in Fornax. Left: Grayscale map (~ 42' per side) of surface 
density of stars with 18.5 <V < 23, from the photometry of [Stetson et al.| ([1998[ reproduced 



with permission from The University of Chicago Press). Notice the crescent- like feature near the 



center. Right: Smoothed stellar map from V, I photometry of Coleman et al. (2005b, reproduced 
by permission of the American Astronomical Society), who detect two shell-like features (one is 
visible ~ 1.3 deg northwest of the center) aligned with lobes along the morphological minor axis. 
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Such peculiarities extend to the ultrafaint satellites as well. Spectroscopic surveys reveal that 
Segue 1 is superimposed on at least one stream of stellar debris ( Geha etal.|2009[ Niederste-Ostholt 



et al. 2009 Simon et al. 201 1 1, offset by ~ 100 km s ^ from Segue 1 in velocity space. Segue 1 also 



shows hints of extended stellar structure superimposed on Sgr debris, although interpretation of 
these extremely low-surface-brightness features in terms of tidal disruption remains controversial 
dBelokurov et al.|2007HGeha et al.|2009t [Niederste-Ostholt et al.|2009| [Simon et al.|201ip . Segue 
2 appears to be embedded in — and comoving with — a stream of stellar debris, perhaps from a 
tidally disrupted parent system CBelokurov et al.[[2009[). Ursa Major II has flattened morphology 



and distorted outer morphology that suggest ongoing tidal disruption ( Zucker et al.|2006at Munoz 



et al. 2010[ ). Another strong candidate for disruption is Bootes III, which has irregular, clumpy 
morphology and a large measured velocity dispersion of a ~ 14 km S-' ( [Carlin et al.|[2009[ ). The 
alignment of stellar distortions in Leo IV and Leo V hints at a low surface brightness 'bridge' 
spanning the ~ 20 kpc between these systems ( [Belokurov et al.|20"08| [Walker et al.|2009a] [de Jong 



et al.[[2010[ note, however, that [Sand et al. [[20 10[ detect no such feature in deep photometry around 



Leo IV). Stars nearest the center of the Willman 1 satellite exhibit near-zero velocity dispersion 
and a mean velocity that is offset from the rest of Willman 1 members by ~ 8 km s"' (Willman 
et al.|20lT] ). Each newly discovered galaxy seems to present a new quirk of its own. 



2.3. Stellar Velocities 



Assuming that dSphs are purely stellar systems truncated by tidal interaction with the Milky 
Way, such that ~ i?K, [Qstriker et al.[ ( [1974[ ) apphed Equation|4]to estimate the mass of the Milky 
Way. Inverting the calculation by assuming an isothermal Galactic halo with Vcirc = 225 km s"\ 



Faber & Lin[ ( |1983| ) estimated masses of dSphs. For the nearest dSphs (Draco, Ursa Minor, Carina, 



Sculptor), Faber & Lin estimated mass-to-light ratios of M/Ly > 10[M/Ly]Q, suggesting dSph 
dark matter. Faber & Lin further used their estimates to predict, via the virial theorem, values of 
> 10 km s~' for the internal stellar velocity dispersions of dSphs. 



At the same time, Aaronson ( |1983[ ) provided the first actual measurement of a dSph's internal 
velocity dispersion. Aaronson used the Multiple Mirror Telescope (MMT, which then consisted of 
six 1.8-meter mirrors working in tandem) to acquire high-resolution {R ~ 30000) spectra for three 
individual carbon stars in Draco. Figure [8] displays the spectra, which were sufficient for Aaronson 
to measure precise line-of-sight velocities of -298.7 ±0.9 km s"' (with a follow-up measurement 
of -297.6 ±0.6 km s"^), -300.2 ±0.6 km s'^, and -279.7 km Aaronson calculated that such 
measurements require, at the 95% confidence level, an intrinsic velocity dispersion of cr > 6.5 
km s"'|^ From dynamical arguments based on the virial theorem ( Illingwor5i|| 1 976 , Richstone & 
Tremaine][1986 . Section 4.1.1 ), such a large dispersion indicates a large dynamical mass-to-light 



^Aaronson added to the final article proof a measurement of -285.6 ± 1.1 km s ' for a fourth, non-carbon star, 
further supporting a large dispersion. 
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ratio M/Ly ^ 30\M /Ly^Q, confirming tlie prediction of Faber & Lin (1983) and indicating the 
presence of dark matter. 

Figure |9] plots the number of stars observed in dSph stellar velocity surveys as a function 
of time. The three decades of observations separate neatly into 'epochs' defined by the available 
instrumentation. 



2.3.1. Small-Number Statistics 



The 1980s yielded the first precise velocity measurements for individual stars in Draco (Aaron- 
son||1983[ ), Carina, Sculptor, and Fornax (Seitze r & Frogel||1985[ ), Leo I and Leo II ( |Suntzeff et al. 



1986| ) and Ursa Minor ( |Aaronson & 01szewski||1987a| ), and follow-up observations of Sculptor 
( [Armandroff & Da Costa|[T986l ) and Draco ( [Aaronson & 01szewski||1987bl ). These samples typ- 



ically included < 10 stars per galaxy and indicated velocity dispersions of a ~ 6 - 10 km s , 
suggesting that dSph dynamical mass-to-light ratios reach at least double the values estimated for 



globular clusters. Strong implications for the particle nature of dark matter (Tremaine & Gunn 



1979, Lin & Faber|1983[ Section 6.2) drew immediate attention. Nevertheless, skepticism regard- 
ing small samples, velocity precision and the unknown contribution of binary orbital motions to 
the measured velocity dispersions demanded further observations and better statistics. 



2.3.2. Confirmation 



In the 1990s, several groups accumulated velocity samples for tens of stars per dSph. Includ- 
ing Aaronson's original work, nine seasons of observations with the MMT eventually produced 
velocity samples that reached ~ 20 members in each of the Draco and Ursa Minor dSphs ( |0T 
szewski et al.|1995 Armandroff et al.|1995 ), yielding velocity dispersion measurements of a ~ 10 



km s and dynamical mass-to-light ratios M/Ly ~ 15\M /Ly^Q for both galaxies. Meanwhile, 
medium-resolution {R ~ 12000) spectra from the William Herschel Telescope gave velocity sam- 



ples for tens of stars in each of Draco, Ursa Minor and the newly discovered Sextans (Hargreaves 
et al.||1994b|a[ |1996bl ). High-resolution spectra from ESO's 3.6-m and NTT telescopes ( |Queloz 
et al. 1995| ) and Keck/HIRES ( Mateo et al.||1998) delivered precise velocities for 23 and 33 stars 
in Sculptor and Leo I, respectively. In all cases, velocity dispersions of a > 6 km s"^ indicated 
M/Ly > 10[M/Lv]q. 



Providing an early demonstration of the efficiency of multi-object fiber spectroscopy, Ar- 
mandroff et al. (1995]) compiled samples of ~ 100 velocities in each of Draco and Ursa Minor 
using the HYDRA multi-fiber spectrograph at the KPNO 4-meter telescope. This data set included 



many repeat measurements, which Olszewski et al. (1996) used to estimate a binary fraction of 
~ 0.2-0.3 for periods of ~ 1 year. Based on Monte Carlo simulations, Olszewski et al.| ( "1996| ) 
and Hargreaves et al.| ( 1996a) concluded that binary motions contribute negligibly to the velocity 
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Fig. 8. — Examples of spectra for individual dSph stars. Left: MMT echelle spectra {R ~ 30000) 
for three carbon stars in Draco, used for the first measurement of a dSph velocity dispersion 
(Aaronson 1983| , reproduced by permission of the American Astronomical Society). Right: 
Keck/DEIMOS spectrum {R ~ 6000) for the brightest red giant in Segue 1 ( |Geha et al.]|2009[ re- 
produced by permission of the American Astronomical Society), with absorption features labeled 
and best-fitting model overplotted. 
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Fig. 9. — Growth of sample sizes available for internal kinematic studies of the Milky Way's 
dSph satellites. Plotted as a function of time is the cumulative number of line-of- sight velocities 
measured for individual stars (including foreground contamination) targeted in dSph spectroscopic 
surveys. References: |Aaronson| ( |1983] ); |Seitzer & Frogell ( |1985[ ); [Suntzeff et al.| ( |1986] ); |Arman 



droff & Da Costa] ( fT986l ); |Aaronson & 01szewski| ( |1987b|aD ; pateo et al.| ( fT99T] ); |Da Costa et al. 



( fT99T] ); |Suntzeff et al.| ( fT993l ); [Mateo eraL] ( [T993] ); |Hargreaves eraL] ( |1994b|a| ); |Armandroff et al. 



( fT995| );|Vogt et al.'(T99?); |Queloz et al.| ( fT99? ); 'Olsz ewski et al.| ( fT997 );'Hargreaves et al. (1996b); 
Mateo et al.| ( |1998| ); |Kleyn a et al.i(l2002','200"3); |Tolstoy et al.|p004l );|Kleyna et al. (2005); ,Munoz] 



eTaL] ( |2006l ); [Munoz et al., (2006) ; ,Battaglia et al.| ( |2006l ); [Westfall et al.. (2006) ; .Walker et al. 



( |2006l |2007aD ; |Koch et al.| p007b|aD ; [Martin et al] ( |2007l ); [Simon & Geha[ ( [2007] ); [Sohn et al. 



(2007); Mateo et al. (2008); Koch et al. (2009); Geha et al. (2009); Walker, Mateo & Olszewski 



(|2009); 'Walker et al.' (2009a); Belokurov et al.'('2009^;'Carli n et aL] ( [2009| );[Willman et al.| ( [20TT] ); 
[Simon et al.|(2011| ); |Aden et al.|(2011) ; |Koposov et al.| ( |2011bl ); [BattaglIa"etal.|(2011] ). 
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dispersions measured for classical dSphs (Section [33]) . 



In the southern hemisphere, Mateo et al. ( 1991 ) used Las Campanas Observatory's 2.5-meter 
telescope and '2D Frutti' echelle photon counter to measure velocities for 44 Fornax stars, in- 
cluding an outer field that showed the same velocity dispersion (cr ~ 10 km s"') as the central 
stars, indicating M/Ly ~ 10[M/Ly]Q. After measuring a velocity dispersion of cr ~ 7 km s"' and 
M/Ly ~ AO\M/Lv]q from the velocities of 17 Carina stars, Mateo et al. ( 1993 ) noted a scaling rela- 
tion among dSphs: the dynamical mass-to-light ratios of classical dSphs are inversely proportional 
to luminosity, suggesting similar dynamical masses of - \QPMq (Section 11). 



2.3.3. Large Samples 



The gap in Figure |9] between 1998-2002 signifies a period not of inactivity but rather of 
construction. During this time, wide-field multi-object spectrographs were built for the world's 
largest telescopes. Over the past decade, surveys with these new instruments have increased stellar 
velocity samples from tens to thousands per dSph. Kleyna et al. (2002 2003} [2004') used the Wide- 
Field Fibre Optic Spectrograph (WYFFOS) at the 2.5-meter Isaac Newton Telescope to measure 
velocities for ~ 100 stars in each of Draco, Ursa Minor and Sextans, respectively. These samples 
were sufficiently large to examine the velocity distribution as a function of distance from the dSph 
center (e.g., [Wilkinson et al.|2002l|Kleyna et al.|2002[[2004HWilkinson et al.|2004[ Section [4|). 



Soon thereafter, the Dwarf Abundances and Radial Velocities Team (DART) used the FLAMES 
fiber spectrograph at the 8.2-meter Very Large Telescope (VLT; UT2) to measure velocities and 
metallicities (derived from the strength of the calcium-triplet absoprtion feature at ~ 8500 A) for 
~ 310, ~ 560 and ~ 175 members of Sculptor, Fornax and Sextans, respectively (Tolstoy et al 



2004 Battaglia et al. 2006 Battaglia et al. 2011). Theses samples yielded the discovery that all 



three of these dSphs contain multiple, chemodynamically independent stellar populations (Figure 



TO[and Section [4X2]). Also with VLT/FLAMES, [Koch et al.[ ( [2007a| ) measured a velocity disper- 
sion of (J ~ 7 km s"^ from ~ 170 members of Leo 11, providing what remains the largest published 
sample for this galaxy. 



2011 



Meanwhile Munoz et al. (2006 ) used archival VLT/FLAMES spectra (see also Fabrizio et al 
to measure velocities for ~ 300 Carina members and added another ~ 45 members from 



spectra obtained sequentially with the MIKE spectrograph at the Magellan/Clay 6.5-meter tele- 
scope. The extra members observed with MIKE extend to ~ 5 times Carina's limiting radius as 
determined from photometry (Figure [6]), indicating that Carina has lost mass to tidal interactions 
with the Milky Way. 



Koch et al. (2007a) used a pair of multi-slit spectrographs — the Gemini Multiobject Spectro- 



graph (CMOS) at the Gemini-North 8-meter telescope and the Deep Imaging Multi-object Spec- 
trograph (DEIMOS) at the Keck 10-meter telescope — to measure velocities for ~ 100 members 
of Leo I. [Sohn et al.[ ([2007]) added another ~ 100 members from their own observations with 
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Fig. 10. — Sculptor's two chemodynamically independent stellar sub-populations (Battaglia et al. 



2008 [ reproduced by permission of the American Astronomical Society ; see also Tolstoy et al. 



2004[ ). Left: Surface brightness profiles for Sculptor's red giants (black), red-horizontal branch 



(red) and blue -horizontal branch stars. Right: Velocity dispersion profiles, calculated for subsets 
of relatively metal-rich (top left) and metal-poor (top right), as determined from the VLT/FLAMES 
spectroscopic sample of Tolstoy et al. (2004). For comparison, the lower- left panel plots the ve- 
locity dispersion profile measured from the composite population. 
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Keck/DEIMOS, and|Mateo et al.|p008[) contributed velocities for ~ 300 Leo I members using the 



Hectochelle multi-fiber spectrograph at the MMT. The latter two studies explored larger radii and 
both found kinematic evidence for tidal streaming motions in the outskirts of Leo I. This result 
is surprising given Leo I's current distance of ~ 250 kpc ( Irwin & Hatzidimitriou|[T995 ), leading 
Mateo et al. ( |2008 ) to suggest that Leo I's orbit is nearly radial. 



Operating from 2004-2011, the Michigan-MIKE Fiber Spectrograph (MMFS; R ~ 20000), 
built by Mario Mateo for the Magellan/Clay 6.5-meter telescope, provided what remain the largest 
homogeneous velocity samples for 'classical' dSphs. The public catalog of Magellan/MMFS ve- 
locities includes ~ 775, ~ 2500, ~ 1365 and ~ 440 members of Carina, Fornax, Sculptor and 
Sextans, respectively ( [Walker, Mateo & Olszewski 2009). Figure 11 displays the Fornax data, 
including sky positions as well as the two quantities measured from each spectrum: line-of-sight 
velocity and a spectral index that indicates the pseudo-equivalent width of the Mg-triplet feature at 
~ 5 170 A. Data from a similar survey conducted in the North with MMT/Hectochelle will soon be- 
come public. Figure [12] displays velocity dispersion profiles calculated from the Magellan/MMFS 
and MMT/Hectochelle data, demonstrating that the luminous regions of classical dSphs have ap- 
proximately constant velocity dispersion. 



2.3.4. (Necessarily) Small Samples 



With kinematic samples for 'classical' dSphs growing exponentially, discoveries of 'ultra- 
faint' dSphs with SDSS data generated a wave of interest in the faintest Milky Way satellites, and 
efforts to obtain spectroscopic follow-up began immediately. Kleyna et al. ( 2005[ ) contributed a 
first result echoing Aaronson's original study of Draco: from Keck/HIRES velocities for 5 mem- 
bers of Ursa Major I, Kleyna et al. (2005) estimate that a > 6.5 km s"^ with 95% confidence. 



Given UMal's low luminosity, simple dynamical models imply M/Ly > 5OO[M/Ly]0. Munoz 



et al. (2006) used the HYDRA multi-fiber spectrograph at the 3.5-meter WIYN telescope to mea- 
sure velocities for 7 members of the Bootes I dSph, obtaining a velocity dispersion of ~ 6.5 km 
s"^ and M/Ly > 13O[M/Lv]0. 

Whereas spectroscopic surveys of 'classical' dSphs target bright red giant branch (RGB) stars, 
the least luminous 'ultrafaint' satellites host few RGBs. Samples for even tens of stars for such 
objects require observations of faint stars near the main sequence turnoff and are feasible only 
with the largest telescopes. [Martin et aL] ( |2007[ ) and [Simon & Geha[ ( [2007 ) used Keck/DEIMOS to 
observe tens of velocities in 10 of the 1 1 'ultrafaints' known at the time, measuring cr > 3 km s"^ 
and concluding that these objects indeed have extremely large dynamical mass-to-light ratios, of 



order M/Ly > 1OO[M/Ly]0 and larger. Geha et al. (2009. see example spectrum in the right panel 



of Figure [8]) and [Simon et al.[ ( [20TT[ ) followed with a Keck/DEIMOS survey of Segue 1, measuring 
(T ~ 4 km s"^ and concluding that this object is the 'darkest galaxy', with M/Ly ~ 34OO[M/Ly]0. 



Aden et aL[ ( [2009[ ) and |Koposov et al.[ ( [201 lb[ ) used VLT/FLAMES to measure velocity disper- 



sions for Hercules and Bootes I, respectively, that while indicative of large dynamical mass-to-light 
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ratios, are both smaller than previously measured with Keck/DEIMOS. [Aden et al.| ( |2009[ ) obtained 
a smaller velocity dispersion after using Stromgren photometric criteria to remove foreground in- 
terlopers. Koposov et al. ( |2011b !) used a novel observing strategy that included ~ 15 individual 



45-60 minute exposures taken over a month. After measuring velocities for each exposure, Ko- 



posov et aL] ( |2011bp were able to resolve binary orbital motions directly and to exclude stars that 



showed significant velocity variability (Section [33 



2.4. The Smallest Galaxies 



Figure [13] displays two scaling relations defined by the structural and kinematic observations 
discussed above, and lets one compare the properties of objects classified as dSphs directly with 
those of objects classified as globular clusters. Parameters for globular clusters are adopted from 
the catalog of |Harris[ ( [T996| 2010 edition). 



The left-hand panel of Figure [13] plots luminosity against size, characterised by the projected 
halflight radius (Belo kurov et~aL][2007t [Gihnore et ai:][2007| [Martin et al.[[2008| [Sand et al.[[20TT] ). 
Here one can appreciate another conclusion of Shapley's (1938b) regarding the first known dSphs: 
'The Sculptor and Fornax systems might be called greatly expanded giant clusters'. Indeed, while 
the luminosity distributions of dSphs and globular clusters overlap substantially, most dSphs have 
Rh ^ 100 pc while most globular clusters have < 10 pc ( [Gilmore et al. 2007 The region 
between 10 < Rh/pc < 100 is populated only by globular clusters with My < -4 or by dSphs with 

Current kinematic results suggest that this separation in luminosity between the smallest 
dSphs and the largest globular clusters is not merely an artifact of classification, but points to a 



fundamental structural difference. The right-hand panel of Figure 13 plots the product Rna^/ (LyG) 
(dimensionally a mass-to-light ratio) against luminosity. In the region of overlapping size, the less 
luminous objects tend to have larger velocity dispersions, amplifying the separation in luminos- 
ity such the smallest, faintest dSphs have the largest dynamical mass-to-light ratios of any known 
galaxies (Kleyna et al.[[2005[ [Munoz et al.[[2006t [Martin et all[2007t [Simon & Geha|[20n7t [Geha 
et al.[[2009^ ,Simon et al.[[201 1 ),^[ Furthermore, the ultrafaint dSphs extend a relation under which 



M31 hosts several 'extended' globular clusters with halflight radii as large as several tens of pc (e.g., Huxor et al. 



2005 1, but no similar population within the Milky Way has yet been discovered. 



^Some ambiguity regarding the masses of the smallest, faintest dSphs results from the convergence of three relevant 
quantities — the typical velocity measurement error, the measured velocity dispersions, and the potential contribution 



to the measured dispersions from binary orbital motions — on the same value, ~ 3 - 4 km s ' ( jSimon & Geha|2007 



[Martin et al.][2007| [McConnachie & C6te|[20T0t [Simon et al.[|20TT] [Koposov et aL]|2011b| l. For some faint dSphs 
the most compelling evidence for large amounts of dark matter comes from stellar-atmospheric chemistry rather than 
kinematics. The faintest objects classified as dSphs tend to have metallicity dispersions (cr[Fe/H] ?i 0-4 dex 



Geha 



et al.|200^ Norris et al.[[2010 Kirby et ar][201 1 Willman et al.|20lT i that are indicative of prolonged and perhaps 



multiple episodes of star formation, thereby requiring gravitational potentials sufficiently deep to retain interstellar 
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less luminous dSphs have larger M/Ly ( Mateo et al.||1993 , Mateo] 1998 SectionjS]). The disconti- 
nuity in dynamical M/Ly between dSphs and globular clusters seems to mark a boundary between 
objects with dark matter and those without. 



3. Stellar Velocity Dispersion as a Proxy for Mass 

In rotating spiral galaxies, the circular velocity at radius r relates directly to enclosed mass via 
^circ = GM(r)/r. Ordered rotation in dSphs is dynamically negligible (see next section); instead, 
dSphs receive support against gravity primarily from the random motions of their stars. Therefore 
the estimation of dSph masses is a fundamentally statistical enterprise. The simplest statistic that 
characterises dSph stellar dynamics is the dispersion of velocities along the line of sight, a. For 
a relaxed system of characteristic size R, the virial theorem implies oc GM/R. In principle, 
measurements of a dSph's size and velocity dispersion provide a simple estimate of its mass (Sec- 
tion]?]). In practice, one must be aware of effects that can inflate measured values of the velocity 
dispersion above equilibrium values. 



3.1. Rotation 



In the simplest case, solid-body rotation about an axis misaligned with the line of sight will 
induce a gradient, dV /dR, in the line-of-sight velocity distribution. For example. Figure [T4| plots 
mean velocity along the axes for which the observed velocity gradients are maximal in Carina, 
Fornax, Sculptor and Sextans. Even at the radii of the outermost observed stars, any ordered 
motion due to rotation is limited to RxaaxdV /dR<3 km s"^ , negligible compared with the observed 
velocity dispersions of cr ~ 10 km s"^ 



For Carina and Fornax the amplitude and orientation of the observed velocity gradients are 
consistent with a perspective effect induced not by rotation, but rather by these dSphs' systemic 
orbital motions transverse to the line of sight ( Kaplinghat & Strigari||2008 , Walker et al.|2008 ), as 
measured from HST astrometry ( Piatek et al.||2002[ 2003 2007[ ). The observed signal in Sculptor 
cannot be attributed to its measured proper motion ( [Schweitzer et al.|ri995[ [Piatek et"aL]|2006| ), 
and thus Sculptor may have a residual rotational component ( Battaglia et al.|2008 1, albeit one that 
contributes weakly (Vrot/o" ^ 0.5) to the measured velocity dispersion. 



media despite pressures generated by stellar feedback. An adequate discussion of the relationships between dSph 
kinematics and stellar chemistry is beyond the scope of the present work; Tolstoy et al. ( 2009 j ) provide an excellent, 
recent review. 
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3.2. External Tides 



All dSphs considered here orbit within the gravitational potential of the Milky Way and are 
therefore susceptible to external influence from tidal forces. Tides can affect the structure and 
kinematics of a given satellite through a variety of mechanisms, including stripping, shocking, 
'stirring' and various orbital resonances. Depending on the strengths of interactions and timescales 
of subsequent relaxation, tides might inflate observed velocity dispersions and mass estimates 
based thereupon. 

In the most extreme scenarios, tides have been invoked to explain dSph velocity dispersions 
without dark matter (e.g., |Kuhn & Milier|[T989| [Kuhnlp93| |Kroupa|p97t [Fleck & Kuhnl[2003 



Metz & Kroupa|20()7| ). Such explanations are most plausible for dSphs that are closest to the Milky 
Way and exhibit elongated stellar structures, e.g., Ursa Major II ( [Zucker et al.|2006a||Munoz et al. 
2ofol ), Hercules ( [Belokurov eraL||2007t [Coleman et"aL|pn07| [Martin & Jin[[20T0l ), and possibly 



Segue 1 (Niederste-Ostholt et al. 2009). However, generalisation of purely tidal mechanisms to 
explain the apparent dark matter content of the entire dSph population does not account without 
contrivance for the wide distribution of dSph distances (~ 30 < Z)/kpc < 250) or for a monotonic 
metallicity-luminosity relation ( [Mateo|1998[[Kirby et al.|2008[|2011| ). 



While tidal stripping involves the transfer of mass from the satellite to the parent outside a 
particular boundar)]^ tidal 'shocking' involves impulsive injection of energy into the satellite as 
it plunges through the disk and/or near the center of the parent system (e.g., Gnedin et al.[[1999 ). 
While the former process removes mass preferentially from the satellite's outer regions, the latter 
process tends to decrease its central density ( Read et aL][2006ap . Early numerical simulations by 
Piatek & Pry or ([1995) and 'Oh, Lin & Aarseth (1995]) suggest that even strong tidal interactions 
do not significantly inflate a satellite's central velocity dispersion, which therefore can remain a 
reliable indicator of dynamical mass. 

Many recent N-body simulations examine specific phenomenology associated with tidal in- 



teractions. For example. Read et al. (2006a) use simulations to demonstrate that the projection of 
tidal streaming motions along most viewing angles tends to cause velocity dispersion profiles to 
increase at large radius, and then argue that the lack of such upturns in the classical dSphs (with 
the possible exception of Draco — see Figure [72]) limits the severity of current disruption events. 



Penarrubia et al. (2008 ) use simulations to show that repeated tidal encounters cause monotonic de- 
clines in equilibrium values (as evaluated at apocenter) of the satellite's central surface brightness, 
velocity dispersion and scale radius. Further simulations by [Penarrubia et aL (20091 indicate that 
when a dark matter halo is present, tides do not generate a clear truncation in the surface brightness 
profile of the bound remnant, as otherwise is prescribed by the self-consistent dynamical model 



^Equation 



by Read et al. 



4 gives the tidal radius for the idealised case of point-mass potentials. Calculations and simulations 
( |2006b ) demonstrate that stars are actually lost from various depths depending on the internal mass 
distributions of satelUte and parent, as well as on the properties of the stellar orbits themselves (e.g., prograde versus 
retrograde with respect to the orbit of the dSph about the Galaxy). 
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of King (1966). Rather, as the bound remnant relaxes after a pericentric encounter, tidal debris 
generates an 'excess' of stars at radii where the local crossing time (Rc/cr ~ 10 Myr at the core 
radius of a typical dSph) exceeds the time elapsed since the encounter. 

Such simulations provide a context for evaluating the tidal origin of breaks and bumps in 
observed surface brightness profiles of individual dSphs (Section 2.2.3 and Figure |2]), and for 



gauging the severity with which tides influence the observed kinematics. In general, one can 
expect tides to have stronger and more enduring influence on the outer regions of dSphs, and for the 
degree of influence on any particular satellite to scale with orbital parameters. A typical dSph with 
Rc ~ 100 pc and a ~ 10 km s"' requires ~ 20 core crossing times to travel a distance of ~ 100 kpc 
(a typical Galactocentric distance) at speed ~ 250 km s"'; the assumption of dynamic equilibrium 
should therefore hold reasonably well out to several core radii for the majority of dSphs. In practice 
one must use all available information about a satelhte's orbit and outer structure to evaluate the 
likely contribution of tides on a case-by-case basis. 



3.3. Binary Stars 



Unresolved binary orbital motions might contribute significantly to the observed velocity 
dispersions of the least massive galaxies. Olszewski et al. (1995]) use 112 independent velocity 
measurements for 42 stars in Draco and Ursa Minor to identify seven stars that exhibit velocity 
variability. Elimination of these stars from their samples has negligible impact on the measured 



velocity dispersions of Draco and Ursa Minor. Using the larger KPNO/HYDRA sample of (Ar- 
mandroff et al. 1995| ) (373 independent velocity measurements for 185 stars), Olszewski et al. 



( 1996) perform Monte Carlo simulations to estimate that the binary frequency for Draco and Ursa 
Minor stars with periods of ~ 1 year is ~ 0.2-0.3 per decade of period. Even though this fraction 
is larger than the one found in the solar neighborhood (e.g., Duquen noy & Mayor||1991 ), further 
simulations by [Olszewski et al.| ( |1996| ) and |Hargreaves et al.| ( |1996a| ) demonstrate that the scatter 
introduced by binaries is small compared to the measured dispersions of a ~ 10 km s"\ and thus 
that binaries do not significantly inflate dynamical masses of the 'classical' dSphs. 

More recently, measurements of velocity dispersions as small as a ~ 3 km s"' for several ul- 
trafaint satelhtes ( [Simon & Gehal[2007t [Martin et al.|[2007t [Walker et al.[[2009at [Aden et al.[[2009 



Koposov et al.[[20lTbl ) have renewed concerns over the possible contribution of binary motions. 
While simulations by Minor et al. ( 2010j ) find that binaries have little effect on dispersions mea- 
sured for systems with intrinsic velocity dispersions cr > 4 km s"\ complementary simulations by 



McConnachie & Cote (2010) demonstrate that for systems with intrinsic dispersions near zero (as 



would be expected for many of the ultrafaints if they contain no dark matter), binaries can inflate 
measured dispersions to values as high as a ~ 4 km s"^ 

It is therefore necessary to verify the extreme mass-to-light ratios of the faintest (and coldest) 
dSphs with repeat spectroscopic measurements that constrain the velocity variabihty of individual 



stars. Simon et al. (201 1 ) present second-epoch Keck/DEIMOS velocity measurements for several 
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of the Segue 1 stars first measured by |Geha et al.| ( |2009l ). They find that one of Segue I's six red 
giants shows significant velocity variability, as do two fainter stars. Including parametric binary 



orbital distribution functions in a Bayesian analysis of Segue I's velocity dispersion, [Simon et al. 

implying a 



(2011, see also Martinez et al.|2011 1 estimate a velocity dispersion of 3.7;^{ | km s ^ 



dynamical mass-to-light ratio of M/Ly ~ 3AQ0{M/Lv]q and reinforcing their previous conclusion 
that Segue 1 is the 'darkest' galaxy known. 



In a separate study, Koposov et al. ( |2011b ) use ~ 15 VLT/FLAMES observations obtained 
over one month in order to resolve binary motions directly among members of Bootes I. Figure 15 
displays independent velocity measurements for two stars as a function of time. While velocities 
for the star in the left-hand panel are consistent with a constant velocity, velocities for the star in 
the right-hand panel change systematically by ~ 10 km s"'. After discarding probable binaries, 
Koposov et al. ( 201 lb| ) find that most members of Bootes I belong to a cold population with 
dispersion a = 2.4;^ °;^ km s~^ (Fig ure 16 ), significantly smaller than previous single-epoch estimates 
of cr ~ 6.5 km s"^ (Munoz et al.||20(j6; Martin et al. 2007). At present the contribution of binaries 



to the velocity dispersions measured for the coldest dSphs is poorly understood, and more multi- 
epoch studies are required. 



4. dSph Masses 



For a coUisionless stellar system in dynamic equilibrium, the gravitational potential, $, relates 
to the phase-space distribution of stellar tracer^ /(r, v, t), via the coUisionless Boltzmann Equation 
(Equation 4- 13c of [Binney & Tremaine|2008| ): 



df df 

^+v-V/-V$-^ = 0. (5) 



Current instrumentation resolves the internal distributions of neither distance nor proper mo- 
tions for dSph stars. The structural and kinematic observations described in section [2] provide 
information only about the projections of phase space distributions along lines of sight, limiting 
knowledge about / and hence also about $. Therefore all efforts to translate existing data sets 
into constraints on $ involve simplifying assumptions. Along with dynamic equilibrium, com- 
mon assumptions include spherical symmetry and particular functional forms for the distribution 
function and/or the gravitational potential. The most useful analyses identify the least restrictive 
assumptions that are appropriate for a given data set. Modern structural/kinematic data contain 
information that is sufficient to place reasonably robust constraints not only on the amount of dSph 
mass, but in some cases also on its spatial distribution. 



^The distribution function is defined such that f{F,v,t)d^xd^v specifies the number of stars inside the volume of 
phase space iPxiPv centered on {x, v) at time t. 



4.1. Amount 



4.1.1. 'Mass Follows Light ' Models 

A common method for analysing dSph kinematics employs the following assumptions: 

1. dynamic equilibrium; 

2. spherical symmetry; 

3. isotropy of the velocity distribution, such that (v^) = (v^) = (v^); 

4. a single stellar component; 

5. the mass density profile, p(r), is proportional to the luminous density profile, u{r) (i.e., M/L is 
constant, or 'mass follows light'). 

Historically, assumption (5) has been adopted when the velocity dispersion profile o{R) is un- 
available. Examples include early analyses of classical dSphs ( Mateo|1998[ and references therein) 
and initial analyses of ultrafaint dSphs (e.g., |Kleyna et al.||2005 , Munoz et"aL]|2006[ [Martin et"ar 



2007, Simon & Geha 20071. Under this assumption, the steeply falling outer surface brightness 
profiles of dSphs (Section 2.2.1 and Figure [2]) motivate the use of dynamical models that allow for 
truncation by external tides (e.g., Michie 1963 ; King||1966 ). Consider, for example, the model of 
King ( 1966[ see also Chapter 4 of Binney & Tremaine 2008), in which the distribution function 
depends only on energy: 

f{e) = Klixvlr'l^ (exp[£/v2] - 1 ) . (6) 

The relative potential, = $o - $, is defined such that £ = \E' - ^v^ > at radii r < R^. Given 
the assumption that mass follows light, this model is fully specified by dimensionless parameter 
\E'(0)/v^ (or equivalently, Rj^/Rc), the core radius R^ and one of either the central velocity dispersion 

(To or central mass density po, which are all related by R^ = ^^J^^ . 

lUingworth ( |1976| ) shows that under assumptions (l)-(5) and a distribution function of the 



form specified by Equation [6} the total mass is given by 

Mtot 



1677y 













km^s"^ 



(7) 



where parameter rj is determined by concentration c = logiQ[i?K/^c]- Following Mateo (1998), 
many authors adopt r/ ^ 8, which is appropriate for the low concentrations characteristic of classi- 
cal dSphs but is not well-constrained for the faintest dSphs. 



More generally, jRichstone & Tremaine (1986]) show that for 'almost any' spherical, isotropic 
system with constant mass-to-light ratio and centrally cored luminosity profile, the dynamical 
mass-to-light ratio is approximately 
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where the half-brightness radius is defined by S(i?hb) = ^^(0) and is often similar (within ~ 25%) 
to 7?,.. Figure [TT] plots dynamical mass-to-light ratios calculated from Equation [8] against dSph lu 



minosity. Masses obtained for the Milky Way's eight classical dSphs are all Mtot ~ IO^Mq (Mateo 



1998). For the ultrafaints, masses range from 10^ < M^x/Mq < 10'^ (Martin et al.[2007, Simon & 



Geha 2007 p . Dynamical mass-to-light ratios increase monotonically with decreasing luminosity, 



ranging from 10 < M/Lv/{M/Lv]q < 1000. 



4.1.2. Does Mass Follow Light? 



When available, empirical velocity dispersion profiles provide simple tests of the assump- 
tions listed above. Dashed lines in Figure \T2\ display best-fitting Kmg](T966} Equation |6]) models 
constructed under assumptions (l)-(5), fit to the surface brightness profiles of [Irwin & Hatzidim^ 
( |1995[ Figure [2]), and normalised to fit the central velocity dispersions. For all of the 



itriou 



Milky Way's classical dSphs, the central velocity dispersions imply large mass-to-light ratios 



M/Ly > 10[M/Lv]q, as plotted in Figure 17 However, the mass-foUows-light models underpre- 



dict the observed velocity dispersions at large radii, which tend to remain approximately constant 
out to the outermost measured points. 

These discrepancies between empirical velocity dispersion profiles and mass -follows -light 
models in Figure [12] imply that at least one of assumptions (l)-(5) is invalid. In fact all are invalid 
at some level. However, various studies indicate that it is unlikely any one of assumptions (l)-(4) 
alone can be the problem. For example: 



1. Simulations by Read et al. (2006a) suggest that tidal disruption (a violation of the equilibrium 
assumption) is more likely to generate rising velocity dispersion profiles rather than the flat profiles 
observed for real dSphs. 



2. Axisymmetric models of Fornax considered by Jardel & Gebhardt (201 H see Section 4.2.2) 
favor dark matter halos that, while violating the assumption of spherical symmetry, extend well 
beyond the luminous component and therefore also invalidate the mass -follows -light assumption. 



3. [Evans, An & Walker (2009() derive analytically an expression for the anisotropy profile, /3a(r) = 
1 - (vg) / (v^) in terms of surface brightness and velocity dispersion profiles. If mass follows light, 
then the flat empirical velocity dispersion profiles tend to imply unphysical values /3a > I. 

4. While recent observations indicate that some dSphs contain at least two distinct stellar sub- 
populations (Sections 2.2.4 and [4.2.2 ), scale radii of the individual sub-populations are sufficiently 
well constrained (e.g., [Battaglia et al.|2006[|Battaglia et al.pOOSj ) that superpositions of two mass- 
foUows-light models continue to underpredict the observed velocity dispersions at large radii. In- 



deed |B^agIiaetaL] ([2008 ) find that the flat velocity dispersion profile they measure for Sculptor's 
more spatially extended subpopulation continues to imply an even more extended dark matter halo. 



Models that allow for sufficiently extended dark matter halos (violating assumption (5) while 
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retaining assumptions (l)-(4)) can provide good fits to surface brightness and velocity dispersion 
profiles simultaneously (e.g., Pry or & Kormendy|1990[ Wilkinson et al.|2()02 ). On these grounds, 
the empirical profiles shown in Figures [2] and 12 provide the strongest available evidence that 
dSphs have dominant dark matter halos that extend beyond luminous regions. 

Some scenarios for dSph formation and evolution — particularly the tidal stirring mechanism 
of Mayer et al. ( 2001a|b , Section 3.2) and the tidal disruption simulations of Munoz et al. ( 2008) — 
tend to produce configurations in which mass approximately follows light. The results discussed 



above seem to rule out this configuration. However, Lokas (2009 1 finds reasonable agreement with 
mass-foUows-light models in Carina, Fornax, Sculptor and Sextans after trimming velocity sam- 



ples in order to remove member stars classified by an iterative mass estimator ( Klimentowski et al 



2007 1 ) as unbound. This result rests in part on a circular argument, as the adopted mass estima- 
tor (jHeisler et al. 1985[ ) is based on the virial theorem, which itself assumes mass follows light. 
However, the same charge of circularity can be brought against the standard kinematic analysis, 
in which the inclusion of stars at large radius in the kinematic analysis implicitly assumes they 
are bound by a sufficiently extended dark matter halo. Thus conclusions regarding the extended 
structure of dSph dark matter halos are generally sensitive to the assumptions employed when de- 
termining which stars to consider or reject in kinematic analyses. More secure is the conclusion 
that dark matter dominates dSph potentials: even mass-foUows-light models require central mass- 



to-light ratios M/Ly > \Q\M/Lv\q in order to fit the central velocity dispersions of dSphs ( [Munoz 
etal.|2008t[Lokas|2009i and Figure [12]). 



4.1.3. Jeans Analysis 

The methods for mass estimation described in the previous section either employ directly 
or are derived from specific distribution functions f(r,v) that correspond to physical dynamical 
models restricted by particular assumptions. Integration of Equation[5]over velocity space provides 
an alternative starting point in the form of the Jeans equations (see |Binney & Tremaine|2008j ). With 
spherical symmetry one obtains 



u dr r r 



(.(v?)) + 2^=-^, (9) 



where z/(r), (v^)(r), and /3a(r) = 1 - (vg)/ (v^) describe the 3-dimensional density, radial velocity 
dispersion, and orbital anisotropy, respectively, of the (stellar) tracer component. The mass profile, 
M(r), includes contributions from any dark matter halo. While there is no requirement that mass 
follow light, there is also no guarantee that a given solution to Equation [9] — even one that fits the 
data — corresponds to a physical dynamical model (i.e., one for which /(r, v) is non- negative). 



Equation [9] has general solution ( [van der Marel|1994 Mamon & Lokas|2005 ) 



GM{s) 
fir) J, ' s^ 



^(^l) = lkl f{s>{s)^^ds, (10) 
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where /(r) = 2/(ri)exp J^.^ (3a(s)s~^ds. Projecting along the line of sight, the mass profile relates 
to observable profiles, the projected stellar density, E(i?) (Figure |2]), and velocity dispersion, a(R) 
(Figure 12), according to ( Binney & Tremaine|2008 ) 



a\R)J:iR) = 2 



dr. 



(11) 



Equation 11 forms the basis for many methods of mass estimation, including parametric (e.g., 
Strigari et al.]|2006l [2008at |Strigari|[20T0t |Koch et al.l|2007b|a| [Battaglia et al.|[20U8[ [Walker et al. 
2007b[|2009bt [Martinez et al.|20"TT] ) and nonparametric (e.g. [Wang et al.pOOS] ) techniques as well 
as algebraic inversion (e.g., Wilkinson et al.|2004|[Gilmore et al.|2007 ). 



All methods based on Equation 1 1 are limited fundamentally by a degeneracy between the 
function of interest, M(r), and the anisotropy profile, (3a{r), which is poorly constrained by ve- 
locity data confined to the line of sightj^ Consideration of a common parametric method helps to 
illustrate this limitation. For example, it is common to assume that the the gravitational potential 
is dominated everywhere by a dark matter halo with mass density profile 



Pir) = Ps 



1 + 



(12) 



i.e., the generalisation by Zhao ( |1996| ) of the Hernquistj ( 1990[ ) profile. Equation 12 provides a 
flexible halo model in the form of a split power-law, with free parameter a controlling the transition 
from index -7 at small radii (r ^ r^) to a value of -(3 at large radii (r ^ r^). 



From spherical symmetry, the density profile specifies the mass profile via 



M(r) = 47r / s p{s)ds 
Jo 

and the surface brightness profile specifies the (deprojected) stellar density profile via 

1 /•°° dR 
u(r) = — 



(13) 
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dR^R^ 



(14) 



Given values for free parameters p^, r^, a, (3, 7 and an assumption about the otherwise uncon- 
strained anisotropy profile]^ Equation [IT] specifies a product S(i?)cr^(i?) that can then be compared 
with observations. 



* Lokas et al. 



[([2005') develop a Jeans analysis that uses higher-order velocity moments (e.g., (v"*)) in order to reduce 
degeneracy between anisotropy (assumed to be constant) and total mass (effectively normalising a cusped mass profile 
that is assumed to have 7 = 1 in the notation of Equation 12 1. 

^Typical assumptions about anisotropy range in simplicity from = or /3a =constant to /3fl(r) = (/3oo-/3o)'"^/('"^ + 
r^) + /?o (e.g., Strigari 2010 1, introducing as many as three new free parameters. 
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Using Sculptor as an example, Figure [T8] demonstrates the degeneracy that is inherent in the 
standard Jeans analysis. Overplotted on Sculptor's empirical velocity dispersion profile are best- 
fitting models obtained under specific assumptions either about the inner slope of the mass-density 
profile (7 = or 7 = 1), or about the amount of velocity anisotropy ((3a = -0.5, /3a = or /3„ = -1-0.3). 
Despite corresponding to a wide range of total masses and mass distributions (bottom panels of 
Figure [T8|), all of these models can provide equivalent fits to the structural and kinematic data. 

While the Jeans analysis provides useful constraints on none of the halo parameters included 



in Equation 12 it does imply a model-independent constraint on a basic dynamical quantity. All 
models that fit the velocity dispersion profile give approximately the same value for the enclosed 



mass near the halflight radius (middle panel of Figure 18), indicating that this quantity is well 
determined by the available data. Indeed, formal confidence intervals derived from Markov-Chain 
Monte Carlo scans of the full parameter space (e.g., |Strigari et al.||2007a| |Strigari|[2010[ [Walker 
et al.|2009b[|Wolf et al.|2010| [Martinez et al.|201 1| ) show a characteristic 'pinch' near the halflight 



radius (Figure 18). 



Equivalent constraints on the enclosed mass at the halflight radius can be obtained more sim- 
ply by solving Equation |9] under the assumptions that (3a = and a(R) =constant. For a Plummer 
surface brightness profile (Equation [3]), one obtains ( [Walker et al.|2009b ) 

M(R,,)=^-. (15) 



Wolf et al. (2010) show analytically that for various surface brightness and anisotropy profiles, the 
tightest constraint on the mass profile (provided that the velocity dispersion profile is sufficiently 
flat) can be approximated by 

M(rs)=^, (16) 

where ri, is the radius at which dlnu/dlnr = -3. For most commonly- adopted surface bright- 
ness profiles, this radius is close to the deprojected halflight radius (i.e., the radius of the sphere 
containing half of the stars), which typically exceeds the projected halflight radius by a factor of 
~ 4/3. Insofar as the assumption of flat velocity dispersion profiles (and the usual assumptions of 
dynamic equilibrium and spherical symmetry) holds, these simple mass estimators can be applied 
even to the relatively sparse kinematic data available for ultrafaint satellites. Walker et al.| ( |2009b| 



see erratum for updated values) and Wolf et al. (2010) tabulate masses obtained from Equations 15 
and [161 respectively, for ~ 2 dozen Local Group dSphs. 



The equivalence of such seemingly crude estimates to constraints from Jeans/MCMC explo- 
rations of parameter space follows from a combination of facts: 1) Equation |9] deals only with 
velocity moments of the phase space distribution function and not with the distribution function 
itself; 2) confinement of empirical velocity distributions to the line-of-sight component yields little 
information about (3a\ 3) therefore the flat velocity dispersion profiles of dSphs effectively reduce 
the available kinematic information to just two numbers, a and the scale radius that characterises 
the adopted surface brightness profile. The information extracted from the Jeans analysis naturally 
amounts to a simple combination of these two numbers. 
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4.2. Distribution 

Some cosmological and particle physics models make specific predictions about how dark 
matter is distributed within individual halos (Section [6]). Insofar as dSphs represent the struc- 
tures most dominated by dark matter and least affected by the presence of baryons, their internal 



dynamical properties provide the most straightforward tests of such predictions. Section 4.1.3 
demonstrates that so long as dSphs have flat velocity dispersion profiles, the standard Jeans analy- 
sis constrains only one number, characterising the amount but not the distribution of dark matter. 
However, the restrictive assumptions employed in the standard Jeans analysis overlook structure 



that is present in the data available for many dSphs (Section [2.2.4[ ). Recent analyses that are devised 
to exploit more fully the available empirical information have been able to place useful constraints 
not only on the amount, but also on the distribution of dark matter within individual dSphs. 



4.2.1. Indirect Constraints 



The presence of stellar substructure within dSphs can provide extra leverage not easily ex- 
ploited in the context of equilibrium dynamical models. For example, it has long been known that 
Ursa Minor's stellar component is 'lumpy' ( iHodge||1964b[ [Olszewski & Aaronson||1985[ |Irwin & 
Hatzidimitriou|1995[|Palma et al.|2003| Section 2.2.A) . While Ursa Minor has a velocity dispersion 



Kleyna et al. (2003) find that a secondary peak in Ursa Minor's stellar density 



of (J ~ 10 km s~' 

field, offest by ~ 20' (~ 400 pc) from the nominal center, exhibits a cold dispersion of < 1 km 
s"'. Interpreting this feature as a loosely bound star cluster captured by Ursa Minor, Kleyna et al. 
( |2003 ) perform N-body simulations to examine how its stability depends on the external gravi- 
tational potential (assumed to be dominated by dark matter) of Ursa Minor. Whereas simulated 
clusters remain intact for a Hubble time when the host potential has a central 'core' with constant 
density (7 = in the notation of Equation [72]) on scales larger than the orbital radius, they disrupt 
on timescales of < 1 Gyr when the host potential has a centrally divergent, or 'cusped' density pro- 
file (7 > |). Therefore, if the observed clump indeed corresponds to a star cluster free of its own 
dark matter component, then its survival provides indirect evidence that the dark matter density of 
Ursa Minor is constant over the central few-hundred pc. 

Substructure in the form of Fornax's five globular clusters provides another example of indi- 
rect evidence for a cored dSph potential. Four of these clusters are projected near (within a factor of 
~ 2) Fornax's halflight radius {Rh ~ 670 pc; Irwin & Hatzidimitriou|1995 ). Hernandez & Gilmore 



( 1998) show analytically that the rate at which the orbits of such clusters decay due to dynamical 



friction depends on the underlying dSph potential. Subsequent numerical simulations by [Sanchez- 
Salcedo et al. ( |2006| ) and Goerdt et al.| (2006) demonstrate that in a cusped potential, dynamical 
friction would require only a few Gyr to bring the clusters from their present positions (assuming 
the projected distances from Fornax's center are not much smaller than the true distances) all the 
way to Fornax's center. On the other hand, in a cored potential dynamical friction would bring the 
clusters only as close as the core radius, where the harmonic potential inhibits further decay. Thus 
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these considerations provide indirect evidence for a core of constant density over the central few- 
hundred pc in Fornax. Further simulations by Goerdt et al. (2010) and Cole et al. ( 2011[ ) suggest 
that the transfer of angular momentum from a sinking cluster to the central dark matter is capable 
of transforming an originally cusped into a cored potential, a possibility that might be relevant for 
the interpretation of evidence for cored dark matter halos in some dSphs (Section [6!!]) . 



4. 2. 2. Constraints from Models 



Several groups have recently developed dynamical and/or kinematic models that exploit struc- 
ture that is present in dSph data but is not considered in the Jeans analyses discussed in Section 



4.1.3[ For example, the discoveries of two chemo-dynamicaUy independent stellar sub-populations 
in Sculptor ( [Tolstoy etal.|2004[ Section|2T4]), Fornax ( [Battagha etal.|2006| ) and Sextans ( |Battaglia 
et al. 2011| ) enable analyses of two tracer components in the same potential. After imposing a 



metallicity cutoff to separate Sculptor's two sub-populations, Battaglia et al. (2008 ) find that cored 
rather than cusped potentials provide better simultaneous fits in a Jeans analysis (Section 4.1.3[ ) of 
the two sets of surface brightness and velocity dispersion profiles. Using the same empirical pro- 
files for Sculptor, I Amorisco & Evans ( |2012| ) confirm this result by modelling both sub-populations 
with anisotropic King-Michie ( King|1962 , Michie|1963[ King|1966 ) distribution functions. These 
dynamical models again favor cored (7 = 0) rather than cusped (7=1) potentials, with a likelihood 
ratio sufficient to reject the hypothesis of cusped potentials with confidence > 99%. 



Jardel & Gebhardt ( |2011 l) take a different approach, constructing axisymmetric three-integral 
Schwarzschild (1979| models for both cored and cusped potentials that also allow for a central 
black hole. For a given potential, libraries of stellar orbits are calculated and each orbit receives 
a weight based on fits to the observed distribution of velocities within discretely binned radii. 
Notice that while the Jeans analysis is sensitive only to the variance of the velocity distribution 
in a given bin, the Schwarzschild method is sensitive to the shape of the distribution. IJardel & 



Gebhardt] (201 1 1 find that their models constructed from cored potentials fit Fornax's velocity data 
significantly better than those constructed from cusped potentials. Within the context of their 
adopted models, they also place an upper limit of < 3.2 x IO'^Mq on the mass of any central black 
hole. Breddels et al. ( 2012| ) independently develop Schwarzschild models for Sculptor, concluding 
that the available data rule out steep cusps (7 > 1.5) but are consistent with slopes in the range 

0<7<lH 



Chaname et al. 



(2008|) have formulated a Schwarzschild method that operates on discrete velocity measurements, 
avoiding the binning procedure altogether Efforts to apply this and other discrete Schwarzschild methods to dSph 
data are underway. As these and other methods for analysing the shapes of dSph velocity distributions continue to 
develop, it will be beneficial and perhaps necessary to allow for contributions from multiple stellar populations. Failure 
to consider distinct populations can generate spurious conclusions regarding orbital structure. To see why this is so, 
consider that the superposition of two Gaussian distributions need not be Gaussian. 
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4.2.3. Direct Measurement 



IWalker & Penarrubia (2011) introduce a method for measuring the slopes of dSph mass 
profiles directly from spectroscopic data and without adopting a dark matter halo model. For a 
given dSph, [Walker & Penarrubia (2011) use measurements of stellar positions, velocities and 
spectral indices to estimate halflight radii and velocity dispersions for as many as two chemo- 
dynamically independent stellar sub-populations. Detections of two distinct sub-populations with 
different sizes provide mass estimates M{Rh) oc RhO^ (e-g-, Equation 15 ) at two different radii in the 
same mass profile, immediately specifying a slope. For Fornax and Sculptor, this method yields 
slopes of r = AlogM/Alogr = 2.61;^o37 and F = 2.955^q^39, respectively, on scales defined by val- 
ues ~ 0.2 < i?/,/kpc < 1 estimated for the halflight radii. These slopes are consistent with cored 
(7 = 0) potentials, for which F < 3 at all radii, but incompatible with cusped (7 > 1) potentials, for 



which F < 2 (Figure 19). 



5. A Common dSph Mass? 



Mateo et al. (1993]) noticed that dSph dynamical mass-to-light ratios increase monotonically 



with decreasing luminosity (left panel of Figure 17). Based on mass-foUows-light analyses (Sec- 
tion 4.1.1 ), the Milky Way's eight classical dSphs trace a relationship Xogy^iM /Ly /{M /Ly^o) ~ 
2.5-1- \ QP Ly /Ly^Q (Mateo 1998). This relation implies that if their stellar mass-to-light ratios are 
M/Ly ~ 2.5 [M/Ly]0, then each dSph is embedded in a dark matter halo of mass ~ IQPMq. Al- 
lowing for extended dark matter halos, |Mateo et al.| ( |l993| ) and [Mateo| ( |1998| ) interpret this value 
as the minimum mass that is associated empirically with dark matter. 

Kinematic studies of the newfound ultrafaint satellites alter this picture slightly. The ultrafaint 
dSphs extend the luminosity floor from My ~ -9 (e.g., Draco, Ursa Minor) to My ~ -2 (e.g.. Segue 
1), nearly three orders of magnitude in luminosity. Applying mass-foUows-light models to their 
velocity data for the least luminous dSphs, [Martin et al.] ( |2007| ) and [Simon & Geha[ ( [2007[ ) estimate 
dynamical masses as small as ~ IO^Mq. While these masses imply extreme mass-to-light ratios 
M/Ly > 1OO[M/Ly]0 that extend the monotonic increase in M/Ly with decreasing luminosity, 
they suggest that this relationship becomes flatter toward the smallest luminosities (right panel of 
Figure [TT]). 



Strigari et al. (2008a) extend to ultrafaint satellites the notion of a common dSph mass by 



considering the amount of mass enclosed within a fixed radius of 300 pc (implicitly assuming 
that all dSphs occupy dark matter halos at least this large). Using a Jeans analysis similar to that 
described in Section [4T3| [Strigari etaLl ( [2008a| ) estimate M300 = M{r < 300pc) ~ WMq for dSphs 
spanning five orders of magnitude in luminosity (Figure[20], left panel). Walker et al.| ( 2009b) take 
a different approach, using the model-independent estimates of M{Ri,) provided by Equation [T5] to 
evaluate the hypothesis that all dSphs are embedded in identical dark matter halos characterised 
by a 'universal' mass profile (Figure [20l right panel). While the scatter of M{Rh) values about a 



single power law, M(r) oc r^, is larger (by a factor of ~ 2) than the scatter expected to arise from 
observational errors, it is similar to the scatter about a common value of M300. 



While it has long been known that dSphs do not fit naturally onto scaling relations defined by 
larger elliptical galaxies (e.g., |Kormendy||1985[ ), the roughly exponential decline of dSph surface 
brightness profiles (Section 2.2. 1 ) hints at a structural relationship to spirals ( [Faber & Lin|1 983';'Lin 
& Faber|1983j , perhaps via an evolutionary mechanism such as tidal stirring (Sections |2. 2. l[,Mayer 



et al. 



2001a|b ). Converting estimates of MiRh) to circular velocities v^jj.^. = GM{Rh)/Rh, 



Walker 



et al. (2010) find that the Milky Way's dsph satellites lie on an extrapolation of the mean rotation 
curve estimated for dark matter in spiral galaxies (Figure 21, left panel), logiQ[ycirc/(kms"^)] 



1.5 + 0.51ogio[r/kpc] ( |McGaugh et~aLl|2007| ). [Donato et al.| ( |2009| ) and |Salucci et al.| ( |2011|) find 



similar results: assuming cored dark matter halos of the form p(r) = p^rlir+rQ) ^(r'^+rg) MBurkert 



1995 1, fits to dSph velocity dispersion profiles extend a simple scaling relationship that is common 
to spiral galaxies ( Kormendy & Freeman 2004 ) and is characterised by logjo[po / (^0pc~^)ro / pc] ~ 
2, independently of luminosity (Figure 21 right panel). 



These scaling relationships have provoked a variety of disparate interpretations. For example. 



Boyarsky et al. (2010a I show analytically that 'secondary infall' of dark matter onto isolated halos 



would cause surface density to vary slowly with halo mass, explaining the apparent constancy of 
PqTq over the range spanned by spirals. They argue that the apparent extension to dSphs agrees with 
results from cosmological N-body simulations (e.g., Springel et al.|[2008) , which produce satellite 



halos with approximately the same surface density. On the other hand, McGaugh & Wolf (2010) 
show that dSphs deviate from another spiral scaling relation — the baryonic TuUy-Fisher relation 
('RTF', |McGaugh et al. 2000) — systematically according to the amount by which their halflight 
radii exceed their tidal radii, provided that the latter radius is estimated using Milgrom's (1983) 
Modified Newtonian Dynamics (a framework that naturally implies the RTF relation). 



In any case, scatter among the Milky Way's dSphs increased when Aden et al. (2009) and 



Koposov et al. ( 201 lb[ ) measured smaller velocity dispersions for Hercules and Rootes I, respec- 
tively, than had previously been estimated (Section [23] ). Furthermore, early kinematic results for 
the dSph satellites of M31 suggest systematically smaller masses (at a given halflight radius) than 
estimated for the Milky Way dSphs ( Kalirai et al.|2010t[Colhns et al.|2010[|201 \) . However, more 
recent, larger data sets suggest that the M3 1 dSphs are more similar to the MW dSphs than previ- 
ously suggested ( ToUerud et al.|[20TT ). Pending more detailed spectroscopic surveys, particularly 
of the least luminous and most distant dSphs, the interpretation of scaling relationships that con- 
nect galaxies across such wide ranges of luminosity and morphological type should proceed with 
caution. 



6. Implications for Cosmology and Particle Physics 

Observations indicate that the Milky Way's known dSph satellites have masses 10^ < M(i?/,) /Mq 
\Qp and that dark matter dominates their internal kinematics at all radii, M/Ly > \0\M / Ly^Q (Sec- 



-35- 



tion 



4.1[). Current observations also provide direct and/or indirect constraints on the internal distri- 



butions of dark matter in three dSphs: Fornax, Sculptor and Ursa Minor (Section 4.2 1. For all three, 
the available evidence indicates central 'cores' of constant density on scales of a few hundred pep] 
Taken at face value, these basic results have implications for broader areas of physics. 



6.1. Cosmology 



Observations of structure on large scales (e.g., as inferred from redshift surveys and anisotropy 
of the cosmic microwave background radiation) seem to require a significant contribution to the 
mass budget from non-baryonic dark matter, ^I^m ~ 0.22 (e.g, Bennett et al.||2003 ; Spergel et al. 



2003 1. The 'cold dark matter' (CDM) cosmological paradigm is built on the hypothesis that the 



dark matter consists of fundamental particles that act like a coUisionless gas after decoupling from 
radiation at non-relativistic speeds shortly after the Big Bang. Small cross sections and low thermal 
velocities allow CDM structure to form and survive at high densities in small volumes, thereby 
enabling the growth of structure on small scales in the early universe. 

Calculations of the matter power spectrum associated with popular 'weakly interacting mas- 
sive particle' (WIMP) candidates for the dark matter (e.g., neutralinos with mass > 10 GeV) 
indicate P{k) oc at small scales until coUisional damping and free streaming finally cause an 
exponential decline on sub-parsec (co-moving) scales ( Green et al.||2004 , Diemand et al.|[2005a ). 
The corresponding halo mass function would be approximately dN(M)/dM oc M"" with a ~ 1.9, 
and a galaxy like the Milky Way would host roughly ~ 10^^ satellites in the form of individual, 
self-bound dark matter 'subhalos', 'sub-subhalos', and 'microhalos' with masses > 10~^Mq 
Pofmann et al.pOOltlDTemand et al.|2005at|Springerer al. 2008). 



Observational requirements derived from the current census of Milky Way satellites seem 
rather modest in this context. A viable dark matter particle needs to accommodate the formation 
and survival of only a few tens (or hundreds when correcting for incompleteness of sky surveys, 
Koposov et al.||2008[[Tollerud et aLl|2008| ) of dark matter halos with masses M > IO^Mq (Section 



|4]) around the Milky Way. These constraints allow for significantly less massive, 'warmer' particle 
candidates (e.g., the sterile neutrino, [Dodelson & Widrow|1994y see recent reviews and discussion 



by, e.g., |Abazajian & Koushiappas|2006||Boyanovsky|2008HBoyarsky et al.|2009b|[Kusenko|2009| ) 
whose longer free-streaming lengths might naturally truncate the matter power spectrum at scales 
more similar to those that characterise the smallest galaxies ( [Gilmore et al.|2007| |Bode et al.|2001 



'in a recent preprint, 



Wolf & Bullock 



2012 



argue that the situation is more ambiguous, and specifically that 



the Sculptor data of Walker, Mateo & Olszewski (20091 favor a dark matter halo that is centrally cusped (7=1) 
rather than cored (7 = 0). Their Figure 1 demonstrates that this conclusion — and more generally the ability of their 
Jeans analysis to yield an apparent constraint on 7 — depends critically on their assumption that the data sample a 
single stellar population with isotropic {(3^ = in the notation of Equation 1 1 1 velocity distribution. This assumption is 



incompatible with the multi-component structure that is present in the spectroscopic data ( Tolstoy et al. 2004,,Battaglia| 
et al.|2006[ [Walker & Penarrubia|20TT] see Sections [ZT4| and p3]i . 
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Maccio & Fontanot|2010||Polisensky & Ricotti|20111|Lovell et al.|2012| ) 



Thus dark matter particle candidates and associated cosmologies can be classified in practicial 
terms according to whether the particles' free streaming plays a significant role in galaxy formation 



(e.g., Boehm et al. 2001 ). For sufficiently massive and 'cold' particles it does not, and other physical 
processes must be invoked to explain the suppression and/or truncation of galaxy formation in 
low-mass halos (e.g., |Klypin et al.|[T999| [Koposov et aL][2009t |Li et al.|[20T0t [Maccio et al.][20T0 



Kravtsov 2010, Font et aL]|2011[ ). The negligible thermal velocities invoked for 'standard' CDM 
particles also imply that N-body simulations can track the growth of structure accurately with 
relatively few particles, making CDM cosmological simulations the simplest, fastest and most 
widely practiced kind. 

Cosmological simulations demonstrate that if gravitational interactions between standard CDM 
particles dominate the formation and evolution of galactic structure, then galaxies ought to be em- 
bedded in dark matter halos that have central cusps characterised by lim^^o p{r) oc r"^, with 7 > 1 
(e.g., IDubinski & Carlberg|[T99T| [Navarro, Frenk & White|[T996| [lW7| [Moore et al.|[T998| [Klypin 



et al.|2001|[Diemand et al.|2005b|[Springel et al.|2008[ ). Observations indicate that most individual 
galaxies with suitable measurements are not embedded in such halos. Instead, rotation curves of 
spiral and low surface brightness galaxies tend to favor dark matter halos with resolved 'cores' 
(7 ~ 0) of constant density (e.g., [Moore|1994t[Flores & Primackll TO.'de Blok & McGaugh|T997 



Salucci & Buricert[[2000t [McGaugh et"aL][200T| [Simon et al.[|2005; Kuzio de Naray et al.[[2006 



2008 [ [de Blok[|2010[ and references therein). These results imply that (standard) gravitational in- 
teractions between CDM particles do not always dominate the formation and evolution of galactic 
structure. 

Indeed galaxies contain baryons prone to interact via forces other than gravity. Many hy- 
drodynamical simulations demonstrate that various poorly-understood baryon-physical mecha- 



et al.[[1996t |E1-Zant et al. 


2001 


Gnedin et al. 


2004, 


Tonini et al. 


2006 


, Romano-Diaz et al. 


2009 


Del Popolo 2010; Govemato et al. 2010; [Pontzen & Govemato 


2011 


; Governato et al. 



2012[ ). Insofar as their baryons are dynamically negligible, dSphs and low surface brightness 
galaxies enable the most direct comparisons to structures formed in CDM-only simulations. In 
this context the available evidence against cusped dark matter halos in Fornax, Sculptor and 



Ursa Minor (Section 4.2) becomes particularly relevant: the viability of standard CDM now re- 
quires that baryon-driven mechanisms can have reduced the central dark matter densities in these 
galaxies to po ^ 5 x lO^A/okpc"^ while leaving behind stellar populations with low luminosities 
10^ < Ly/Li/,0 < 10' and central surface brightnesses 23 < po/(mag/arcsec^) < 25 

Recent work identifies several mechanisms that might accomplish this feat on dSph scales 



Boylan-Kolchin et al. 



(2011) identify a similar (perhaps the same) structural problem, noting that the most massive 
'subhalos' produced in the Aquarius CDM simulation ( Springel et al.[[2008 i have central densities larger than those 
estimated for any of the known dSphs. 
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by invoking either the dynamical coupling of the dark matter to energetic baryonic outflows (e.g., 
Read & Gilmorel[2005| [Mashchenko et al.||2006l [2008t |de Souza et"aL][20TT] ) or the transfer of 
energy/angular momentum to dark matter from massive infalling objects (e.g., 'S anchez- S alcedo] 
et al.||2006[ |Goerdt et aL]|20061 [20T0t |Cole et aL]|201 Hydrodynamical simulations by .Sawala] 



et al.| ( |2010[ ) and [Parry et aL] ( |2011[ ) indicate that the former category of solutions has difficulty 
reproducing other dSph observables — specifically, star formation histories as well as luminosity 
functions and metallicity distributions. The latter category of solutions is difficult to evaluate ob- 
servationally, as the evidence can literally be destroyed (e.g., by tidal disruption); furthermore it 
seems unlikely that such infall mechanisms generate cores of sufficient size{^ Alternatively, cos- 
mological simulations that consider 'warmer' particle candidates demonstrate that the associated 
suppression of small-scale power can naturally (i.e., without invoking baryon physics) produce 
halos with large cores; however, such scenarios seem to require fine tuning of the relative contri- 
butions from various production mechanisms in order to reproduce simultaneously the number of 
observed MW dSphs ( [Polisensky & Ricotti|20TT||Macci6 et al.|2012a|bl ). 



In any case, the emerging challenge for the standard CDM paradigm is not that empirical 
evidence against cusped dark matter halos necessarily rules out the hypothesis that CDM particles 
constitute the dark matter. The poorly understood complexities of baryon physics — along with the 



freedom to invoke other processes, e.g., self-scattering of CDM particles ( [Spergel & Steinhardt 
lOOOHLoeb & Weiner|[20TT||Vogelsberger et al.||20T2l )— leave sufficient flexibility for CDM to be 
rendered consistent with virtually any realistic observation of galactic structure. In fact that is the 
problem. CDM escapes falsification of perhaps its most famous prediction only by withdrawing the 
prediction. While this circumstance does not imply that CDM is incorrect, it does mean that CDM 
currently fails to make accurate predictions regarding the stellar dynamics of galaxies, a primary 
piece of evidence for dark matter in the first place. In this context a decisive outcome favorable 
to standard CDM seems to require the detection of either 1) gravitational interactions involving 
dark matter halos on sub-galactic scales (e.g., via microlensing or perturbations of loosely bound 
luminous structure) or 2) nongravitational interactions involving cold dark matter particles. 



6.2. Particle Physics 



It has long been recognised that the small sizes and large mass densities of dSphs place strong 
constraints on the particle nature of dark matter. For example, Liouville's theorem requires that 
the phase space densities of light, neutral lepton species do not increase after decoupling from 



radiation in the early Universe. Tremaine & Gunn ( 1979) point out that this constraint, combined 
with the necessity that il^ < ^rnMer, places a conservative upper limit on the neutrino mass that is 



^For example, 



Goerdt et al. 



( 2010 1 conclude that a sinking object of mass induces core formation inside a radius 
In this scenario the sinking of Fornax's five surviving globular clusters 

550pc) ~ 5 X IO^Mq and M(~ 900pc) ~ 



where the enclosed halo mass is M(rcoie) ^ M. 

(Ms ^ IO^ Mq) cannot have formed the core inferred from estimates M( 
2 X IO'^Mq (IWalker & Penarrubia|201 1|. 
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summarily violated by lower limits from phase space densities inferred for galaxy halos. Therefore 
neutrinos are not the dark matter in galaxies. This exclusion is most evident on small scales, where 
small volumes demand heavy particles in order to satisfy phase-space requirements. 

For example, using Aaronson's (1983) initial measurement of Draco's velocity dispersion. 



Lin & Faber ( 1983) derive a lower limit of niy > 500 eV. Lake ( 1989) points out that this constraint 



is sensitive to the dubious assumption that mass follows light (Section 4.1.1 1. Gerhard & Spergel 



( 1992) strengthen the argument by turning it around, noting that for more viable neutrino masses 
of ~ 30 eV, the core radii of dSph halos would need to be unrea listically large (> 10 kpc) 
to accommodate model-independent lower limits of po ^ O.OSMqPC"^ ( Pry or & Kormen(^ |l990) 



on their central densities. Generalising the phase-space argument oflTremaine & Gunn| ( |1979 ) to 
relativistically decoupled warm dark matter candidates, Dalcanton & Hogan,(2001] ) show that of all 



galaxies, dSphs provide the most stringent limits, > 700 eV and > 300 eV for thermal and 



degenerate fermions, respectively. Most recently and more specifically, Boyarsky et al. (2009a) 



use phase-space arguments to conclude that > 1.7 keV if the dark matter consists of sterile 
neutrinos produced via non-resonant mixing with active neutrinos. 

Any positive identification of a dark matter particle will require the detection of its non- 
gravitational interactions. Experiments at the Large Hadron Collider might find evidence for such 
interactions, as might the various experiments designed to detect directly the scattering of dark 
matter particles in Earth's orbital path. Alternatively, high-energy photons might be released if 
dark matter self-annihilates ( |Gunn et al.]|1978| |Stecker||1978[ ) or decays ( |Pal & Wolfenstein||1982 



Boyarsky et al.|2006 ; Kusenko|2006| ), providing an opportunity for indirect detection. 



Their large mass-to-light ratios, low astrophysical backgrounds and close proximities make 
the Milky Way's dSph satellites popular targets in the search for annihilation and/or decay products 
(e.g., [Evans et al.|2004||Strigari et al.|2008bt[Kuhlen|2010| ). For annihilation, the differential 7-ray 



flux (units cm ^ s ' sr ^ GeV ^) received on Earth in solid angle AQ is given by 



_ 1 (crv) dN^f 



X y(A^]), 



(17) 



where is the particle mass, (crv) is the (velocity-averaged) cross section, dN^/dE^ is the energy 
spectrum of products and 



7(Afi) = 



An 



p^{l,VL)dldVL. 



(18) 



This '7-factor' represents the astrophysical contribution to the signal and is specified by the integral 
of the squared dark matter density, p^(Z, Vt), over line of sight / and solid angle Vt. The equation for 
the flux due to decay events is similar, except that the integral is taken over the dark matter density 
raised only to the first power. Published constraints on J come directly from constraints on p{r) 



obtained in parametric Jeans analyses of the sort described in Section 4.1.3 and demonstrated in 
Figure [18] (e.g. , [Strigari et al.|2007b| [Martinez et al.|2009t [Charbonnier et al.|2011| ) . 



At present, dSph surveys conducted with atmospheric Cherenkov telescopes (e.g., Pieri et al. 


2009 Essig et al.|2009 ; H. E. S. S. CoUaboration et al.|201 1 


the VERITAS collaboration: Vivier et 
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aL][20TT| [Aleksic et aL][20TT] ), x-ray (e.g., |Boyarsky et al.|[2007t |Loewenstein et al.|[2009| |Riemer 



S0rensen & Hansenl|2009[ [Loewenstein & Kusenko|[20T0| [Boyarsky et al.|[20T0bt |Loewenstein & 



Kusenkol |2012[ ) and gamma-ray telescopes (e.g., |Abdo et al.||2010[ [Scott et aL]|2010^ [Ackermann 



et al.,,2011 1 yield no unambiguous detections!^ From Equation 17 upper limits on photon flux 



translate into upper limits on the cross section (crv) for a given particle mass and annihilation 



channel. For example, Figure 22 plots 95% upper limits on (crv) derived from Fermi-LAT ob- 
servations of Milky Way dSphs, based on two years of data fro m the planned five-year mission 
( Geringer-Sameth & Koushiappas|201 1 ; Ackermann et al.'201 1). Dotted lines at (av) ~ 3 x 10"^^ 



crtrs 



mark the 'generic' cross section expected for WIMPs with mass ~ 0.1 - 1 TeV (e.g.. 



Jungman et al.|1996 Feng|2010 1. WIMPs having this combination of mass and cross section would 
have decoupled from radiation with relic abundance r2wiMP ~ 0.2, the value cosmology requires of 
the dark matter (a coincidence sometimes referred to as the 'WIMP miracle'). For particle masses 
below < 100 GeV, the combination of kinematic and high-energy data available for dSphs is 
now encroaching upon the cross section most readily associated with WIMPs{^ Over the next 
decade, searches for dark matter and/or its byproducts will intensify with large-scale efforts at 
existing facilities and with new instrumentation that will provide unprecedented sensitivity (e.g., 
CTA Consortium|2010D . 



The status of dark matter as a particle will depend critically on the outcomes of direct and 
indirect dark matter detection experiments that are either ongoing or planned for the near future. 
Also at stake is the motivation and context for studying dark matter phenomenology in dSphs. An 
unambiguous, positive detection of dark matter byproducts emitted from dSphs would establish 
the existence of a new particle and would provide a unique means for measuring its mass and cross 
section; via Equation [17} constraints on such parameters would be only as good as dynamical 
constraints on the dark matter density profile. In case of detections in other objects, e.g., the 
Galactic center (see ' Hooper & Linden||201 1[ ) or galaxy clusters (see Han et al.|2012 ), high-energy 
and dynamical constraints from dSphs would provide important consistency checks in the regime 
of minimal astrophysical background. In the case of direct detection and characterisation of a new 



^' Loewenstein & Kusenko 



(2010 1 interpret a Chandra detection of monochromatic 2.5 keV) emission from the 



direction of the Willman 1 satelhte as a possible signal of sterile neutrino decay. However, Boyarsky et al. (2010b I 
argue that non-detections of this feature in the Galactic halo, M31 and several other dSphs rule out a dark matter 
origin. Indeed, [Loewenstein & Kusenko] ( [201 2| l report no detection of the 2.5 keV feature in follow-up XMM- 
Newton observations of Willman 1; corresponding limits on the mass/mixing angle of sterile neutrinos depend on how 
reh ably the 'irregular' stellar kinematics of Willman 1 ( Willman et al.|20lT 



Charbonnier et al. 



(2011 



Section 



2.2.4 1 trace its mass. 



I use published kinematic data to estimate less stringent limits of (crv) < 10"^^ cm^s"' (at 
~ 10 GeV, cf. Figure|22]) for individual dSphs. Possible reasons for this discrepancy include different assumptions 
about the dark matter halo profile ( Geringer-Sameth & Koushiappas|20lT and Ackermann et al.|201 1 1 adopt J values 
previously estimated under the assumption that dSph dark matter halos follow NFW (7 = 1 in the notation of Equa- 
tion [T2]i profiles; [Charbonnier et al.|201 1| estimate J values by marginalising over uncertain halo shape parameters), 
different assumptions about the energy spectrum ( |Geringer-Same th & Koushiappas 20 lT| and [Ackermann et al.|2011| 
explicitly consider annihilation via bb and t'^t~ mechanisms; Charbonnier et al.[201 1 [ consider a conservative spectrum 
averaged over a variety of plausible annihilation channels) and/or different assumptions about detector sensitivity. 
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particle in the laboratory, dSph phenomenology would help to establish or rule out association with 
cosmological dark matter. In case of failure to detect non- gravitational interactions involving dark 
matter particles, consideration of the extreme phenomenology exhibited by dSphs would help to 
inform alternative explanations for dark matter. 



7. Some Considerations for Future Work 

Topics and/or questions for future investigations include the following: 

• The least luminous galaxies are often disproportionately interesting. What sorts of bound 
stellar/gaseous objects will the next generation of sky surveys detect? 

• This field will continue to develop in proportion to 1) the rate at which data sets grow in 
size and quality and 2) the efficiency with which analytical tools can extract the information 
contained therein. Promising strategies for analysis include the formulation of statistical 
likelihood functions that depend on discrete measurements (e.g., of individual stellar po- 
sitions and velocities) rather than on binned profiles, and exploitation of multi-component 
structure/substructure that is present in the data. 

• All dSphs clearly deviate from the simple models invoked to characterise their observed 
properties. In order to understand dSphs as galaxies, one must consider not only the struc- 
tural and kinematic data emphasised above, but also chemical abundances, star formation 
histories (i.e., stellar masses), internal substructure and external environment. 

• The outer stellar structures of dSphs carry valuable information about the gravitational com- 
petition between the Milky Way and its satellites. Combined with deep photometric surveys 
that map the Galactic stellar halo, wide-field spectroscopic surveys that reach large galacto- 
centric distances in dSphs will reveal details of the transition from satellite to host potential. 

• The standard CDM hypothesis seems to require baryon-physical processes to have formed 
large constant-density cores in galactic dark matter halos. Do the internal mass distributions 
of galaxies correlate with other observables directly related to baryon physics (e.g., lumi- 
nosity, metallicity, star formation history)? Constraints on core sizes for the least luminous 
dSphs would be extremely valuable. 

• Can observations and/or experiments provide clear evidence for non-gravitational interac- 
tions involving a dark matter particle? Can observations provide evidence for gravitational 
interactions involving dark matter halos on subgalactic scales? 

• Are there cosmological and/or particle physics models for dark matter from which can be 
extracted accurate predictions regarding the stellar dynamics of the most dark-matter domi- 
nated galaxies? 
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Fig. 11. — Magellan/MMFS spectroscopic data for Fornax (Walker, Mateo & Olszewski 2009 



updated to include the full sample of ~ 3200 members). Left: Sky positions of individual stars. 
Right: velocities and spectral indices (pseudo-equivalent widths of the Mg-triplet absorption fea- 
ture. Color indicates membership probability, as estimated from position, velocity and spectral 
index distributions. The ellipse indicates the limiting radius, Rk (Equation [T]). 
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Fig. 12. — Velocity dispersion profiles observed for the Milky Way's eight 'classical' dSphs 
HWalker et al.|200 7aVMateo et al. |2008H Walker et al.|2009bD . See also 'Kleyn a etaL] ( |2002l [2004| ) ; 
[Wilkinson et al.| ( ,2004) ; ,Munoz'eraL] p005, ,2006] ); |Sohn et al.| ( |2007) ; ,Koch et al.| ( |2007b|aD ; 



Battaglia et"aL] p008[ |2011| ). Overplotted are mass-foUows-light | King| ( |1966| ) models (Section 



4.1. 1| ) normalised to reproduce the observed central dispersions. Failure of these models to repro- 



duce the large velocity dispersions at large radius provides the strongest available evidence that 
dSphs have dominant and extended dark matter halos. 
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Fig. 13. — Left: Luminosity versus size (updated from [Belokurov et al.||2007| [Gilmore et al. 



2007[ [Martin et al.]2008j ), for Milky Way satellites including objects classified as globular clusters 



(open triangles; data from |Harris| ( |1996[ 2010 edition)) and dwarf spheroidals (filled circles with 
error bars; data from compilations by [Irwin & Hatzidimitriou||1995| |Mateo|[T998| [Martin et al. 



2008, Sand et al. 2011). The apparent lack of objects toward low luminosity and large size is a 



selection effect that reflects the surface-brightness limit of the SDSS survey ( [Koposov et al.|2008] ). 
Right: Dynamical mass-to-light ratio (modulo a constant scale factor) versus luminosity. Updated 
from |Mateo eraL] ( [T993l ); |Mateo| ( fT998] ); [Simon & Geha| ( |2007] ); |Geha et al.| ( |2009l ). Dotted lines 
correspond to constant masses of lO^M©, lO^M© and \QPMq. 
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Fig. 14. — Line-of-sight velocity gradients in the Milky Way's 'classical' dSph satellites. Panels 
display mean velocity as a function of distance along the axis corresponding to the maximum ve- 
locity gradient, from the samples of Walker, Mateo & Olszewski] ( 2009] ). The observed gradients 
give rise to maximum amplitudes Rj„axdV / dR < 3 km s"\ significantly less than the velocity dis- 
persions (cr ~ 10 km s"'). Using VLT spectra from the DART survey, |Battaglia et al. (2008) report 
similar results for Sculptor, measuring a gradient of 1 ■(i^22. km s"Meg"' along the morphological 
major axis. 
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Fig. 15. — Direct detection of binary stars from individual velocity measurements in Bootes I 
( [Koposov et al. 2011b, reproduced by permission of the American Astronomical Society). Left: 
foreground star with velocities consistent with a constant value. Right: member star ((v)boo ~ 105 
km s"^) exhibiting a velocity increase of ~ 10 km s"^ over one month. See section 3.3 for a 
discussion of the effects of binary orbital motions on measurements of dSph velocity dispersions. 
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Fig. 16. — Velocity distribution of the Bootes I dSph, from the VLT/FLAMES data of [Koposov 
et al. ( 201 lb[ reproduced by permission of the American Astronomical Society), not including 
1 1 stars observed to have significant velocity variability. In the best-fitting two-Gaussian model, 
70% of the stars belong to a 'cold' component with a ~ 2 km s"^ and the rest belong to a hotter 
component with cr ~ 9 km s"^ 
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Fig. 17. — Dynamical mass-to-light ratio, derived from mass -follows -light models (Section 4. 1 . 1 1, 
versus luminosity. Data for the Milky Way's 'classical' dSphs (filled circles) are from the re- 
view of [Mateo (1998). Data for the 'ultrafaint' dSphs (open squares) are from Keck/DEIMOS 
observations by Simon & Geha ( 2007| ) and Martin et al. ( 2007| ). The dotted line corresponds to 
[M/Lv/[M/Lv]q] = 1.5 + lQp /{L/Lv^q) (Section|5j). 
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Fig. 18. — Degeneracy in the Jeans analysis of dSph mass profiles (Section 4.1.3). In the top 
panel, overplotted on Sculptor's empirical velocity dispersion profile are best-fitting cored (7 = 
in the notation of Equation [12]) and cusped (7=1) density profiles, as well as models that assume 
isotropic (/?„ = 0), radially anisotropic (/?„ = +0.3) and tangentially anisotropic {(3a = -0.5) velocity 
distributions. All of these models provide equivalent fits to the data. The bottom two panels plot 
mass and density profiles corresponding to each model. Dotted curves indicate 95% confidence 
intervals, as determined from a Jeans/MCMC analysis that assumes constant anisotropy and lets 



the unspecified parameters in Equation 12 vary freely (Section 4.1.3). The vertical dotted line 
indicates Sculptor's projected halflight radius. 
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Fig. 19. — Empirical constraints on the slopes of mass profiles in Fornax and Sculptor, based 
on estimates of M(r/,) for each of two chemo-dynamically independent stellar sub-populations in 
each galaxy (Section 4.2.3 Walker & Penarrubia|201 1[ reproduced by permission of the American 
Astronomical Society). Points in the left two panels indicate constraints on the two halflight radii 
and masses enclosed therein, with lines indicating the maximum slopes allowed for cored (7 = in 
the notation of Equation [T2| T = AlogM/Alogr < 3) and cusped (7 = 1, F < 2) overplotted. The 
right-hand panel indicates posterior probability distributions for the slope in each galaxy. The data 
rule out the cuspy (7 > 1, F < 2) profiles that characterise the cold dark matter halos produced 
in cosmological simulations (e.g., Navarro, Frenk & White 1997) with significance > 96% and 
> 99% in Fornax and Sculptor, respectively. 
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Fig. 20. — A common mass? Left: mass enclosed within a radius of r < 300 pc, estimated from 
a Jeans analysis (See Section [4.1.3| ), versus luminosity ( Strigari et al.|2008a| reproduced with per- 
mission). Right: mass enclosed within the projected halflight radius, estimated from Equation [T5j 
versus halflight radius ( Walker et al.|2009b| reproduced by permission of the American Astronom- 
ical Society). Arrows indicate the trajectory followed by a satellites as they lose up to 99% of their 
original stellar mass, from N-body simulations by |Penarrubia et aL (2008 ). Overplotted are various 
dark matter halo models ( [Walker et al.|2009bl ). 
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Fig. 21. — Similarity with spiral galaxies (Section |5])? Left: circular velocities (vdrc = 
GM{Rh) / Rh) for the Milky Way's dSph satellites lie on the extrapolation of the mean rotation 



curve due to dark matter in spiral (McGaugh et al. 2007) and low surface brightness (Kuziode 



Naray et al. [2008 1) galaxies ( Walker et al.|2010l reproduced by permission of the American Astro- 



nomical Society). Right: Dark matter surface densities inferred from fits of Burkert ( 1995) halos 
have the same constant value inferred from spiral galaxy rotation curves (reproduced from Dwarf 
Spheroidal Galaxy Kinematics and Spiral Scaling Laws, by P. Salucci et al., MNRAS, in press, by 
permission of John Wiley & Sons Ltd.). 




Fig. 22. — Exclusion of WIMP self-annihilation cross sections, based on Fermi-LAT non- 
detections (2-year data) of gamma-rays in the Milky Way's dSph satellites (reprinted with per- 
mission from Geringer-Sameth & Koushiappas, Phys. Rev. Lett. 107, 241303 (2011; left) and 
Ackermann et al. (The Fermi Collaboration), Phys. Rev. Lett. 107, 241302 (2011; right). Copy- 
right 2011 by the American Physical Society). 



